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CHAPTER  1  INTRODUCTION 


The  conventional  edge  emitting  laser  diode  uses  either  a  single-stripe 
or  multiple  stripe  structure  consisting  of  an  active  layer  in  a  plane 
parallel  to  the  semiconductor  surface;  two  cleaved  facets,  functioning  as 
laser  cavity  mirrors,  are  perpendicular  to  this  active  layer.  It  is 
relatively  difficult  to  fabricate  and  diagnose  in  situ,  and  only 
one-dimensional  laser  arrays  have  been  monolithically  fabricated.  The 
surface  emitting  laser  (SEL),  on  the  other  hand,  utilizes  both  surfaces  of 
the  wafer  as  laser  mirrors  and  the  light  output  is  perpendicular  to  the 
wafer  plane.  This  makes  planar  fabrication  and  diagnosis  easier,  and 
more  importantly,  allows  for  two-dimensional  (2-D)  monolithic 
integration  and  array  formation.  This  2-D  SEL  array  can  be  scaled  in 
number  and  size  while  maintaining  a  high  filling  factor.  The  SEL  array 
diodes  can  also  be  individually  addressed  and  controlled,  and  the  output 
beam  can  be  easily  accessed  from  both  the  front  and  back  surfaces  of  the 
wafer.  These  features  are  very  important  and  desirable  when  a  large 
array  of  coherent  emission  semiconductor  lasers  is  needed.  ^ 

In  addition,  the  SEL  is  typically  very  short  in  active  length  in 
comparison  to  an  edge  emitting  laser  diode  (about  10  mm  versus  200  to 
300  um  ).  As  a  result,  single  mode  or  single  frequency  operation  can  be 
achieved  relatively  easily  from  SELs.  Because  the  optical  cavity  is  small, 
the  photon  round  trip  time  is  short  (about  0.2  psec),  making  very  high 
frequency  modulation  possible.  Because  of  the  large  emitting  area  of 
SELs  (typically  10  um  diameter,  as  compared  with  1  um  in  conventional 
edge  emitting  laser  diodes),  the  output  beam  divergence  of  a  SEL  is 
expected  to  be  smaller  than  that  of  an  edge  emitting  laser.  These  are  all 
favorable  factors  for  application  in  high  data  rate  optical  communication 
and  signal  processing. 

Although  the  SEL  promises  to  be  very  attractive  for  2-D  arrays  and 
for  various  applications  that  require  fast  response  and  single  frequency 
operation,  the  SEL  device  technology  is  still  relatively  young  and  not 
fully  established.  So  far  most  of  the  SEL  device  development  effort  has 
been  conducted  in  Japan  by  K.  Iga  and  coworkers,  and  few  research 
efforts  had  been  conducted  in  the  United  States  until  recently,  after  our 
research  team  at  LMSC  began  development  work  in  collaboration  with 
the  University  of  California  at  Berkeley. 
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The  SEL  structure  reported  by  the  Japanese  group  utilized  a  vertical 
junction  similar  to  that  in  edge  emitting  laser  diodes,  without  any 
transverse  current  confinement  and  with  relatively  poor  optical 
confinement.  As  a  result,  the  device  exhibited  a  very  high  threshold 
current  (about  a  few  hundred  mA).  In  addition,  the  laser  cavity  is 
formed  by  an  etched  mirror  and  ring  electrode,  which  has  relatively 
poor  optical  quality.  Therefore,  new  and  novel  device  structures  and 
configurations  are  needed  to  improve  the  SEL  performance  both  in 
output  efficiency  and  threshold  current. 

The  realization  of  a  2-D  SEL  array  by  the  Japanese  group  was  limited 
to  a  2  x  2  array  with  only  75  percent  yield,  and  a  larger  SEL  array  has  so 
far  not  been  demonstrated.  In  order  to  prove  the  concept  that  SELs  can 
be  easily  scaled  into  a  2-D  array,  a  research  effort  in  this  area  is  needed. 

The  main  objective  of  this  program,  therefore,  is  to  develop  a  novel 
SEL  design  and  structure  with  improved  performance  characteristics, 
and  to  demonstrate  the  feasibility  of  scaling  the  SEL  to  a  2-D  array.  This 
report  summarizes  accomplishments  and  progress  toward  this  objective 
during  the  previous  18-month  period.  Chapter  2  describes  the  technical 
concepts  and  approaches  we  have  taken  in  designing  and  fabricating  the 
SEL.  Chapter  3  presents  the  results  in  device  design  and  fabrication,  laser 
characterization,  array  development  and  theoretical  analysis.  Chapter  4 
summarizes  the  overall  progress  and  accomplishments  and  discusses  the 
technical  issues  remaining  to  be  solved.  In  the  last  chapter,  we 
recommend  a  follow-on  program  that  allows  these  remaining  technical 
issues  to  be  fully  investigated.  This  effort  will  bring  SEL  technology  to  a 
more  mature  stage,  and  allow  its  true  potential  and  practicality  to  be 
further  explored. 
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CHAPTER  2  TECHNICAL  CONCEPTS  AND  APPROACHES 


2.1  Key  design  considerations 

The  two  prerequisites  for  efficient  operation  of  a  laser  diode  are  good 
optical  and  current  confinement.  These  two  features  have  been 
successfully  implemented  in  the  conventional  edge  emitting  laser  diodes. 
However,  the  SEL  diodes  reported  so  far  do  not  provide  either  type  of 
confinement.  Therefore  a  special  structure  must  be  devised  to 
implement  these  confinement  features.  We  chose  a  transverse  buried 
heterojunction  structure,  in  which  the  current  is  injected  laterally  into 
the  active  region.  This  structure  not  only  improves  the  current  and 
optical  confinement  but  also  provides  a  longer  and  more  uniform  active 
length  to  assure  higher  gain  and  power  output. 

Another  consideration  for  this  new  transverse  junction  SEL  is  the 
optical  cavity  structure.  The  need  for  a  low  loss  optical  cavity  is  quite 
evident  from  a  threshold  gain  calculation.  The  threshold  gain  of  a  laser  is 
given  by 


9th  =  «  +  (1/L)ln(1/R) 

where  a  is  the  internal  losses  due  to  the  free  carrier  absorption  and 
scattering,  L  is  the  active-layer  length,  and  R  is  the  reflectivity  of  the 
laser  cavity  (assume  both  facets  have  equal  reflectance).  For  a 
conventional  edge  emitting  laser  diode  of  L  =  250  um,  R  =  0.32,  and  a  = 
10  cm**,  the  threshold  gain  is  estimated  to  be  about  60  cm"*.  To  keep 
gth  low  in  an  SEL  with  a  short  active  length,  we  require  a  cavity  mirror 

reflectance  of  R  ~  0.95  (cavity  length  L  =  10  mm,  gth  =  60  cm'*,  and 
same  internal  losses). 

To  achieve  this  high  a  mirror  reflectance  in  the  SEL  cavity,  two 
approaches  were  explored.  One  approach  was  to  form  alternate  layers  of 
GaAs  and  AlGaAs  material  with  quarter-wavelength  thickness  within  the 
laser  active  region,  or  on  the  top  and  bottom  portion  of  the  laser  with  the 
active  region  in  between.  Another  approach  was  to  deposit  dielectric 
coatings  on  the  output  end  of  the  SEL  instead  of  the  alternate 
quarter-wavelength  GaAs/AlGaAs  multilayer.  Either  approach  can  in 
theory  form  a  laser  cavity  with  a  reflectivity  as  high  as  99  percent. 
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For  SELs,  it  is  extremely  important  that  the  lasers  can  be  made  of 
uniform  quality  in  layer  thickness,  doping  concentration,  and  cavity 
length.  The  first  two  requirements  can  be  met,  and  the  third 

requirement  can  be  facilitated  by  using  either  the  molecular  beam 
expitaxy  (MBE)  or  metal  organic  chemical  vapor  deposition  (MOCVD) 
growth  technique.  Uniform  cavity  length  is  achieved  by  using  a  selective 
etching  technique  that  stops  at  a  clearly  defined  GaAs/AIGaAs  interface. 
These  approaches  will  assure  formation  of  a  mirror-like  flat  surface 
upon  which  dielectric  coatings  can  be  deposited  to  give  superior  mirror 
quality. 

2.2  Basic  surface  emitting  laser  configurations 

Based  on  the  design  considerations  discussed  previously,  we 

conceived  a  SEL  utilizing  a  transverse  buried  heterojunction  (TBH) 
structure  with  a  vertical  laser  cavity.  Two  basic  SEL  configurations  of 
this  TBH  structure  are  shown  schematically  in  Fig.  1.  One  is  a  distributed 
feedback  (DFB)  type  and  another  is  a  distributed  Bragg  reflector  (DBR) 
type  resonator  configuration. 

The  DFB  type  SEL  has  a  smaller  gain  volume  that  can  operate  at  low 
threshold  current.  It  also  requires  no  dielectric  mirror  and  has  good 
mode  selectivity.  However  this  structure  can  have  nonuniform  spatial 
gain  and  its  performance  is  relatively  sensitive  to  the  DFB  layer 

thickness.  The  DBR  type  SEL  on  the  other  hand  is  easy  to  design,  is  less 
sensitive  to  fluctuations  in  layer  thickness,  and  it  can  be  adjusted  for 

maximum  gain  or  output  by  varying  the  active  region  volume.  It  also 
has  a  more  uniform  spatial  gain  profile.  However,  the  DBR  structure 
requires  a  high  reflectivity  dielectric  mirror  coating  or  multiple  layers  of 
GaAs/AIGaAs,  and  typically  needs  a  high  operating  current. 

We  chose  the  DFB  type  structure  as  the  main  approach  during  this 
project  because  of  the  encouraging  results  obtained  from  this  structure 
earlier,  and  because  of  the  relative  simplicity  in  device  fabrication.  The 
detailed  design  and  fabrication  process  of  this  DFB  type  SEL  and  its  array 
are  described  in  the  following  section. 
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CHAPTER  3  PROGRESS 


3.1.  Device  design  and  fabrication 

The  successful  fabrication  of  any  new  device  usually  requires  utilizing 
well-developed  processing  technology  but  often  involves  combining 
several  processing  steps  in  a  new  sequence.  The  initial  processing  steps 
must  be  compatible  with  and  must  not  be  affected  by  the  processing  to 
follow.  Because  a  semiconductor  laser  requires  the  synergistic  action  of 
electronic  and  optical  structures,  the  demands  on  performance  and  on 
the  nature  of  the  fabrication  steps  is  even  greater.  The  following 
sections  describe  the  complex  series  of  steps  required  to  fabricate  our 
SELs,  discuss  the  issues  involved  in  our  choice  of  electro-optical 
structure,  and  present  the  results  from  our  fabricated  samples. 


3.1.1.  DFB  type  SEL  structure 

As  we  described  previously,  the  main  SEL  structure  we  developed 
during  this  project  is  a  DFB  type  cavity  as  shown  in  Fig.  1.  In  this  case, 
the  performance  of  the  optical  cavity  is  determined  by  the  sequencing  of 
the  epitaxially  grown  multilayer.  The  effective  reflectance  of  the  cavity 
structure  increases  with  the  number  of  pairs  of  layers,  and  with  the 
refractive  index  difference  between  the  two  layers  within  a  pair.  Each  of 
the  layers  is  nominally  one  quarter  wave  of  optical  thickness  at  the 
design  (Bragg)  wavelength.  The  resonator  loss  is  least  at  the  Bragg 
mode,  and  increases  for  wavelengths  greater  or  less  than  this  mode. 
Since  a  symmetric  DFB  structure  with  a  periodic  real  refractive  index 
modulation  will  not  lase  precisely  at  the  Bragg  wavelength,  an 
asymmetry  must  be  introduced  into  the  resonator  to  induce  oscillation  of 
the  Bragg  mode.  The  "phase  shift  layer",  which  in  our  structure  is  an 
additional  quarter-wave  layer  of  GaAs  located  20  pairs  from  the  output 
facet,  serves  to  provide  this  asymmetry.  It  forces  single  frequency 
oscillation  at  the  design  Bragg  wavelength,  and  also  increases  the  ratio  of 
forward  to  backward  emitted  power. 
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We  chose  an  index  guided  waveguide  structure  for  confinement  of  the 
optical  mode.  Index  guided  structures  have  several  advantages  over 
gain  guided  structures  for  laser  diodes,  including  a  lower  threshold 
current,  a  sharper  turn-on  at  threshold,  and  a  more  stable  near  and  far 
field  radiation  pattern.  Most  other  SEL  designs  (with  the  exception  of  in 
plane  surface  emitters  relying  on  second  order  distributed  Bragg 
reflection  to  eject  the  beam  perpendicular  to  the  wafer  plane)  have 
utilized  gain  guiding  for  optical  mode  confinement.  The  buried 
heterojunction  structure  created  by  the  LPE-reerowth  process  serves  to 
provide  this  real  refractive  index  waveguiding. 

Finally,  we  chose  a  transverse  buried  heterostructure  (TBH)  type  of 
pn  junction  for  electron  and  hole  confinement.  This  type  of  charge 
carrier  control  has  been  used  quite  successfully  in  edge  emitting  lasers, 
and  it  is  easy  in  principle  to  implement.  The  biggest  advantage  a  TBH 
type  of  junction  for  SELs  is  its  compatibility  with  other  planar  processing 
technology. 

3.1.2.  Fabrication  procedure 

The  starting  point  for  the  fabrication  of  a  SEL  is  the  epitaxially  grown 
multilayer  structure.  MOCVD  was  used  to  deposit  the  alternating  layers 
of  GaAs  and  A1Q  3GaQ  ?As  on  a  semi-insulating  GaAs  substrate.  The  wafer 

surface  was  then  masked  and  a  wet  chemical  etch  was  applied  to  create 
a  mesa.  This  mesa  forms  the  core  of  the  optical  waveguide  as  well  as  the 
active  region  of  the  laser.  The  wet  etch  typically  results  in  a  tapered 
profile.  The  mesa  as  shown  in  the  scanning  electron  microscope(SEM) 
photo  (Fig.  2)  is  about  4  mm  by  20  mm,  with  the  smaller  dimension 
increasing  to  about  12  mm  at  the  base  . 


Liquid  phase  epitaxy  (LPE)  is  then  used  to  regrow  an  n-type 
AlQ4Ga06As  cladding  layer  around  the  mesa,  completing  the  optical 

waveguide  structure.  The  silicon  nitride  on  the  top  of  the  mesa  acts  as  a 
mask  for  the  regrowth,  allowing  selective  epitaxy  to  occur  only  around 
the  islands.  Another  nitride  mask  is  then  applied,  and  a  lateral  pn 
junction  is  formed  by  zinc  diffusion  through  an  opening  in  the  n-type 
cladding  region.  The  diffusion  is  timed  so  that  the  diffusion  front  stops 
in  the  active  region,  providing  for  carrier  confinement  in  the  waveguide 
core. 
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In  the  final  steps,  ohmic  contacts  are  applied;  Au/Cr  is  used  to 
contact  the  p  side  of  the  junction,  and  Au-Ge-Ni  is  used  for  the  n  contact. 
The  wafer  is  diced  into  chips  roughly  1.5  mm  on  a  side,  and  the 
individual  chips  are  mounted  in  a  chip  carrier  with  either  a  thermal 
epoxy  or  indium  solder.  Finally,  ultrasonically  bonded  aluminum  wire  is 
used  to  connect  the  chip  carrier  pads  to  the  bonding  pads  on  the  laser 
array. 


Figure  2.  SEM  picture  of  the  cross-section  of  a  DFB  type  SEL 
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Two  fabrication  issues  should  be  discussed  at  this  time.  The  proper 
timing  of  the  zinc  diffusion  for  the  p  side  of  the  junction  is  rather  crucial 
to  the  successful  performance  of  the  laser.  The  diffusion  front  must  stop 
in  the  center  of  the  multilayer,  so  that  the  carriers  are  confined  in  the 
center  of  the  waveguide  core  and  maximum  gain  is  provided  to  the 
optical  mode.  The  position  of  the  zinc  front  varies  from  laser  to  laser, 
depending  on  the  local  wafer  temperature  and  the  local  crystal  quality. 
The  poor  rectifying  performance  of  some  of  the  lasers  may  be 
attributable  to  the  diffusion  front  stopping  on  the  interface  between  the 
multilayer  and  the  regrown  AlGaAs  cladding,  increasing  the  likelihood 
that  carriers  will  leak  through  the  junction  via  interfacial  defects. 

A  second  issue  is  the  problem  we  have  encountered  with  the  adhesion 
of  the  ohmic  contact  metallization.  This  is  caused  by  two  competing 
requirements;  the  window  in  the  silicon  nitride  mask  must  be  small  to 
restrict  the  area  of  the  zinc  diffusion,  but  must  be  large  enough  to 
provide  a  sufficient  area  for  adhesion  to  the  GaAs  surface  (the  adhesion 
to  the  nitride  is  poor).  We  expect  to  eliminate  this  problem  with  a 
different  choice  of  contact  metallurgy.  These  two  issues  also  affect  the 
fabrication  yield  of  our  SELs  and  SEL  array,  leaving  only  a  limited 
number  of  SEL  array  samples  available  for  performance  characterization. 

3.2.  Laser  performance  characterization 

3.2.1.  Device  test  apparatus 

The  first  test  that  was  performed  on  the  newly  fabricated  SELs  was  to 
record  the  current  versus  voltage  response  (I-V)  as  a  check  of  the 
quality  of  the  pn  junction  and  of  the  carrier  confinement.  A  Tektronix 
curve  tracer  was  used  for  this  measurement.  Next  the  near  field 
radiation  pattern  and  the  optical  output  power  versus  current  input 
(L-I)  were  recorded.  A  manual  semiconductor  probe  station  was  useful 
for  examining  the  near  field  at  magnifications  to  500x;  we  modified  a 
Vickers  Instruments  microscope  to  allow  probe  testing  of  our  SELs  and 
SEL  array  samples,  visual  inspection  of  the  near  field  pattern,  and 
recording  of  the  L-I  data.  This  microscope  had  an  accessory  that  utilized 
shearing  techniques  to  measure  image  features  with  roughly  1  um 
accuracy  and  was  used  to  make  dimensional  measurements  on  the 
fabricated  samples.  The  near  field  pattern  was  recorded  using  an  RCA 
infrared  vidicon  with  an  image  digitizer. 
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Lastly,  the  spectral  emission  profile  was  measured  for  selected  lasers. 
A  McPherson  monochromator  was  used  in  combination  with  a  cooled 
germanium  pin  detector  (North  Coast  Model  EO-817L).  The  germanium 
detector  was  chosen  for  its  low  noise  equivalent  power  as  well  as 
reasonably  flat  spectral  response  near  900  nm.  A  helium  cryostat  was 
also  used  to  take  temperature-dependent  data. 


3.2.2.  Electrical  characteristics 


For  all  the  devices  fabricated  on  a  wafer,  the  uniformity  in  the  I-V 
response  was  generally  quite  good.  A  typical  example  of  the  DC 
electrical  response  is  shown  in  Fig.  3.  The  turn-on  voltage  is  about 
1.3- 1.5  volts,  and  the  breakdown  voltage  is  about  -7  volts.  This  result 
indicates  that  the  pn  junction  is  formed  properly  in  the  device.  However, 
some  junctions  exhibit  a  large  amount  of  leakage,  and  others  show  no 
rectifying  behavior  at  all.  We  have  also  observed  a  long  term 
degradation  of  the  pn  junction  and  laser  performance;  this  may  be 
caused  by  migration  of  defects  into  the  active  region. 


3.2.3.  Near  and  far  field  radiation  patterns 


The  refractive  index  difference  between  the  active  region  and 
cladding  layer  of  our  SEL  waveguide  structure  is  about  0.12  and 
therefore  should  be  more  than  adequate  for  providing  strong  mode 
confinement.  This  is  exemplified  by  a  typical  well-confined  near  field 
radiation  pattern  as  shown  in  Fig.  4.  All  of  the  lasers  that  we  examined 
have  exhibited  this  strongly  confined  emission  pattern,  which  also 
indicates  that  the  TBH-type  pn  junction  is  providing  good  carrier 
confinement.  The  light  emission  is  restricted  to  the  rectangular  surface 
emitting  area  in  both  directions. 
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Figure  3.  Typical  I-V  reponse  of  the  SEL  diode 


The  far  field  emission  pattern  observed  by  a  video  camera  shows  a 
near  circular  pattern  with  a  divergence  angle  of  about  8  degrees  as 
shown  in  Fig.  5.  However,  some  laser  diodes  show  a  quite  different  far 
field  when  measured  with  a  rotating  slit  method.  The  emission  is  quite 
broad  in  extent,  which  could  be  due  to  the  strong  spontaneous  emission 
component  observed  at  high  operating  currents. 
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3.2.4.  L-I  characteristics 

I 


For  semiconductor  lasers  of  any  type  the  variation  in  optical 
performance  from  one  laser  to  the  next  can  be  rather  large.  This  is 
caused  by  the  differences  in  crystal  quality,  layer  uniformity,  processing 
steps,  and  in  the  quality  of  heat-sinking.  As  a  result  the  SEL  diodes 
exhibited  a  broad  range  of  optical  performance.  We  show  in  Fig.  6  a  set 
of  L-I  curves  for  an  arbitrarily  chosen  group  of  lasers.  We  have 
measured  a  best  case  output  efficiency  just  under  1%,  and  have  observed 
powers  as  high  as  0.5  mW  from  many  emitters.  In  measurements  of  the 
polarization  of  the  output  from  several  lasers,  we  have  observed  only  a 
slight  difference  in  the  relative  power  in  the  two  orthogonal  polarization 
directions.  This  result  is  not  completely  unexpected;  Because  the 
equivalent  active  layer  thickness  of  the  SEL  is  about  5  mm,  which  is  well 
within  the  regime  in  which  TE  and  TM  mode  reflectances  are 
comparable.  It  is  possible,  however,  that  residual  birefringence  could  act 
to  select  one  preferable  polarization. 


Figure  6.  Typical  L-I  curves  of  SELs 
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The  temperature  dependence  of  the  L-I  response  has  also  been 
measured  and  the  result  is  shown  in  Fig.  7.  The  "characteristic 
temperature"  is  a  quantity  that  is  usually  quoted  to  quantify  the 
performance  of  the  laser  at  various  temperatures,  where  it  is  assumed 
that  the  threshold  current  I[h~  exp(T/TQ)  (T  is  the  temperature  and  TQ 

is  the  characteristic  temperature).  From  this  data  we  determine  that  TQ 
=  1 90K. 
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CURRENT  (mA) 


Figure  7.  Temperature  dependence  of  L-I  response 
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3.2.5.  Emission  spectrum 

The  spectral  emission  from  these  lasers  has  several  components  as 
shown  in  Fig.  8.  At  very  low  currents  (typically  just  a  few  mA)  a  narrow 
peak  appears  above  the  background,  centered  at  about  884  nm.  This 
mode  corresponds  to  the  fundamental  Bragg  mode  as  designed,  and  it 
grows  in  relative  intensity  until  much  higher  currents  are  reached 
(50-100  mA),  when  the  modes  at  shorter  wavelength  begin  to  compete 
for  the  gain.  The  linewidth  of  this  mode  is  somewhat  broader  than 
expected,  and  it  is  difficult  if  not  impossible  to  diagnose  the  reason  for 
this  broad  emission  width.  Possible  explanations  include  thickness 
errors  in  the  epitaxially  grown  multilayer  as  well  as  carrier-induced 
shifts  in  the  refractive  index,  which  both  degrade  the  finesse  of  the 
resonator,  thus  broadening  the  emission  linewidth. 


Figure  8.  Typical  SEL  emission  spectrum 
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Additional  modes  appear  at  shorter  and  longer  wavelengths;  the 
theoretical  model  described  below  predicts  the  positions  of  these 

longitudinal  modes  with  reasonable  accuracy.  From  an  examination  of 
the  axial  intensity  profile  of  all  the  modes  it  seems  likely  that  these 
modes  can  contribute  to  the  output  because  they  can  utilize  the  spatial 
gain  that  is  unused  by  the  main  Bragg  mode.  This  envelope  spatial  hole 
burning  effect  is  predicted  to  be  less  severe  in  distributed  Bragg 
reflector  (DBR)  lasers. 

3.3.  Laser  array  development 

3.3.1.  Array  design  and  configuration 

Because  of  the  planar  structure  and  surface  emission  of  SELs, 

formation  of  1-D  and  2-D  arrays  is  relatively  simple  to  accomplish;  a  3  x 
3  SEL  array  might  look  like  that  shown  in  Fig.  9  .  By  proper  design  of 

the  photolithographic  masks,  the  SEL  array  can  be  fabricated  using  the 

same  sequence  of  processing  steps  that  is  used  in  the  fabrication  of 
individual  SELs. 


Figure  9.  Schematic  of  a  3x3  SEL  array  configuration 
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Using  the  DFB  type  SEL  diode  structure  described  previously,  we 
fabricated  a  number  of  arrays  with  1-D  and  2-D  configurations.  The 
array  chip  was  mounted  on  a  standard  multiple  pin  chip  carrier  for 
testing  as  shown  in  Fig.  10. 


Figure  10.  Picture  of  a  SEL  array  mounted  on  a  chip  carrier 


Fig.  11  shows  the  photomicroscopic  picture  of  a  2-D  SEL  array  with  6 
x  16  elements  on  a  GaAs  chip  of  about  2  mm  x  2  mm.  The  light  emitting 
area  is  adjacent  to  the  shiny  rectangular  region,  which  is  a  gold  electrode 
pad.  The  spacing  between  each  SEL  element  is  about  100  mm,  limited 
by  the  electrode  pad  size  (about  80  mm  by  80  mm).  The  pad  size  was 
chosen  for  reasons  of  reliability;  a  smaller  pad  delaminates  from  the 
substrate  during  ultrasonic  wire  bonding.  The  light  emitting  area  varies 
from  3  mm  x  15  mm  to  3  mm  x  20  mm. 
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Figure  11.  Photomicroscopic  view  of  a  6  x  16  SEL  array 


3.3.2.  Array  performance 

Fig.  12  shows  the  emission  pattern  of  the  SEL  array  under  cw 
operation.  Because  of  the  different  emission  areas  of  SEL  elements  in  the 
array  and  possible  nonuniformity  in  the  wafer  fabrication  process,  the 
threshold  current  for  each  SEL  varies  widely  between  a  few  mA  to  a  few 
tens  of  mA.  Therefore,  the  array  emission  pattern  was  taken  by 
adjusting  the  driving  current  of  each  SEL  to  provide  a  nearly  uniform 
emission  intensity. 

As  can  be  seen  from  the  emission  pattern,  the  overall  array  has  a 
fairly  good  yield  with  about  95%  of  the  SELs  in  the  array  emitting 
radiation,  and  only  a  few  SELs  out  of  about  a  hundred  do  not  emit  light. 
The  non-emitting  devices  often  exhibited  no  rectifying  behavior, 
indicating  that  the  lack  of  light  emission  was  due  to  little  or  no  carrier 
confinement.  The  carrier  confinement  is  strong  in  all  lasers  that  do  emit 
light,  as  evidenced  by  the  tight  optical  mode  confinement. 

This  result  is  very  encouraging  for  future  array  development. 
Additional  improvements  in  array  uniformity  will  be  forthcoming  as  the 
wafer  quality  and  fabrication  technology  are  refined. 
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Figure  12.  Emission  pattern  of  a  6  x  16  SEL  array 


3.4.  Theoretical  modeling  and  analysis 

In  conjunction  with  the  experimental  work  an  effort  of  theoretical 
modeling  was  undertaken.  The  objective  was  twofold:  (1)  to  aid  in 
interpreting  the  experimental  data  on  fabricated  SELs,  and  (2)  to  provide 
guidance  in  the  design  of  the  laser  optical  resonator.  An  eigenmode 
theory  based  on  a  Bloch  wave  expansion  of  the  electric  field  was 
developed  and  applied  to  the  analysis  of  DFB  type  SELs.  This  treatment 
was  restricted  to  the  below  threshold  regime.  An  alternative  formalism, 
using  thin  film  discrete  propagation  techniques,  was  also  constructed. 
This  model  is  capable  of  describing  strong  mode  coupling  between  the 
counter-propagating  waves  in  the  resonator  as  well  as  treating  gain 
saturation  and  laser  operation  above  threshold.  These  models  are 
described  in  detail  in  the  sections  that  follow. 
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3.4.1.  Bloch  wave  analysis  of  surface  emitting  lasers 

The  distributed  feedback  laser  can  be  described  quite  simply  if 
one  assumes  that  the  solutions  of  Maxwell's  equations  for  an  infinite 
periodic  structure  approximate  the  solutions  for  the  finite  laser 
structure.  The  laser  electric  field  must  satisfy  the  Helmholtz 
equation 


d2E(z,X)  9 

- [g(z,X)  -  i  p(z,X)]2E(z,X)  =  0 

dz^ 


(1-1) 


where  E(z,\)  is  the  electric  field,  g  is  the  material  gain  coefficient,  and 
(3  is  the  propagation  constant.  If  the  gain  and  refractive  index  are 
assumed  periodic  with  period  A, 

g(z,X)  -  i  P(z,X)  =  gQ(X)  -  i  p0(X) 


Q-oo 

+  2^(gq(l)  -  i  Pq(l))exp(i2q7tz/A)  (1-2) 

q=-oo 

then  the  solutions  for  the  electric  field  are  of  the  Floquet-Bloch  type: 
E(z)  =  A(z)  exp(Tz)  +  B(z)  exp(-Tz)  (1-3) 

where 

r  =  G  -  iK  (1-4) 

G  is  the  "effective  gain"  and  K  is  the  "effective  propagation  constant," 
and  A(z)  and  B(z)  are  periodic  functions  of  z.  A(z)  corresponds  to  a 
forward  propagating  mode  and  B(z)  corresponds  to  a  backward 
propagating  mode.  When  the  condition  for  Bragg  reflection  i  s 
satisfied,  i.e.  when  (3q  is  approximately  (the  Bragg  propagation 

constant),  the  waves  scattered  at  each  interface  between  dissimilar 
media  constructively  interfere  and  the  overall  reflectivity  becomes 
high.  If  we  further  expand  A  and  B  in  a  Fourier  series 
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A(z,X)  =  Aq  +  £  Aq(X)exp(i2qKz/A) 

(1-5) 

q  =  -o° 

q=oo 

B(z,X)  =  B0+  £  Bq(X)exp(i2q7rz/A) 

q=-oo 

(1-6) 

we  can  derive  first  order  expressions  for 

the  amplitude  components 

iK+o 

A+P(X)  ~  G+g0+i(8+5rff)  A0 

(1-7) 

B-PW  =  G+g0+i(8+8eff)  B0 

(1-8) 

where 

8  -  Kb  -  P0(X) 

(1-9) 

8eff  =  K„  -  K 

(1-10) 

K+p  *  P+p  +  ‘S+p 

(1-11) 

K.p  =  P.p  +  ig.p 

(1-12) 

We  also  obtain  a  characteristic  equation 

(G  +  iSeff)2  =  (gQ  +  i8)2  +  K+pK.p 

0-13) 

which  is  used  to  determine  eigenmode  solutions  for  the  pair  of 
quantities  (G,  8eff).  If  we  apply  the  above  to  our  DFB  surface  emitter, 

we  obtain  the  results  shown  in  Fig.  13,  where  G,  8eff,  and  8  are 
plotted  as  a  function  of  wavelength  for  the  case  when  the  material 

gain  coefficient  g  =  2000  cm'*.  It  is  seen  that  the  effective  gain  is  a 
maximum  in  the  region  of  the  Bragg  wavelength,  indicating  that  the 
structure  is  working  properly  to  provide  maximum  feedback  at  the 
Bragg  wavelength. 


< 
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We  have  also  used  this  model  to  predict  the  emission  spectrum  of  the 
DFB  type  SEL.  Briefly,  the  laser  was  modeled  as  a  spontaneous 
emission-driven  resonator;  the  formalism  described  above  was  used  to 
determine  the  amount  of  wavelength-dependent  feedback,  and  the 
spectral  variation  of  the  spontaneous  emission  power  was  included.  As 
shown  in  Fig. 14  the  computed  mode  positions  and  the  shift  of  emission 
wavelength  with  temperature  agree  quite  well  with  the  experimental 
data,  but  the  estimate  of  the  relative  emission  intensities  is  rather  poor. 
This  may  be  due  to  the  linear  nature  of  the  model,  which  neglects  gain 
saturation  and  is  therefore  applicable  only  below  laser  threshold.  For 
analysis  of  laser  emission  properties  above  laser  threshold  we  developed 
an  alternative  model  described  below. 


X  (nm) 

Figure  13.  Wavelength  dependence  of  G  and  8eff 
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3.4.2.  Discrete  propagation  model  for  DFB  and  DBR  lasers 


Coupled  mode  theory  is  the  standard  approach  to  modelling 
distributed  feedback  (DFB)  lasers,  and  it  has  proved  to  be  useful  in 
the  first-order  design  of  a  variety  of  semiconductor  laser  structures. 
It  should  strictly  be  used  to  analyze  edge-emitting  lasers,  since  the 
coupled-mode  propagation  equations  include  an  implicit  assumption 
of  weak  interaction  between  the  forward  and  backward  travelling 
waves  in  the  resonator.  SELs  operate  in  the  strong  coupling  regime; 

typically  the  distributed  feedback  coefficient  K  ~  5000-10000  cm'* 

for  surface  emitters,  and  K  ~  50  cm'*  for  edge  emitters.  Thus  some 
loss  in  accuracy  would  ensue  in  applying  the  methods  of  coupled¬ 
mode  theory  directly  to  SELs.  In  order  to  take  this  strong-coupling 
effect  into  consideration  and  to  include  the  effects  of  gain  saturation, 
we  have  extended  the  discrete  propagation  techniques  of  Bjork  and 

Nilsson*  and  have  used  this  model  to  a  analyze  SELs  with  a  DFB 
vertical  cavity  and  with  a  distributed  Bragg  reflector  (DBR)  cavity. 
The  methods  of  Ref.  1  are  particularly  amenable  to  computer 
calculation  and  well  suited  to  analyzing  resonators  having  a 
multilayer  structure.  They  are  also  flexible  enough  to  be  used  in 
cases  involving  facet  reflections,  phase  shift  layers,  amplifying 
regions,  and  dielectric  coatings. 

The  equations  we  used  to  model  these  lasers  comprise  two  sets. 
The  first  set  treats  the  free  propagation  of  the  forward  and  backward 
waves  within  a  layer  according  to 

Aj(Zj+lj)  =  Aj(Zj)  exp[i(Jjlj  +  (1/2)  gj(Aj,  Bj)lj)  (2-1) 

BjUj+lj)  =  Bjlzj)  expl-ipjlj  -  (1/2)  gj(Aj,  Bj)^)  (2-2) 

gj(A.;.  Bj)  =  g0/(l  +  lAjUj)!2  +  IBj(Zj)l2)  (2-3) 

where  Aj  and  Bj  are  the  normalized  forward  and  backward  wave 
amplitudes  in  the  i-th  layer,  Zj  is  the  axial  position  of  layer  i,  (3j  is  the 
propagation  constant,  lj  is  the  i-th  layer  thickness,  and  gj  and  gg  are 

the  saturated  and  unsaturated  intensity  gain  coefficients 
respectively.  The  fields  are  propagated  across  an  interface  between 
layers  i  and  j  using  the  second  equation  set 
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OUTPUT  POWER  NORMALIZED  TO  SATURATION  POWER 


Figure  14.  Temperature  dependence  of  emission  spectrum 


Figure  15.  Power  output  vs  pumping  parameter 
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(2-4) 


(Pi  +  Pj)  (Pi-Pi) 

aj  -  ¥-2*^  Ai  +_ipT Bi  1 

(Pi  -  Pi)  (Pj  +  Pj) 

Bi  =  Fji  [  ia-  A;  +  - 


J  »J  1  2 pj  “i 


2P 


J 


F.  ,£l(i/2)  - 2ii - 

1J  Pi  (HijHjj)^/2) 


Hjj  =  £  1  j  dxdy  WjWj1 


(2-5) 

(2-6) 


(2-7) 


where  wj  is  the  waveguide  mode  field  function  in  layer  i,  and  the 
sum  in  Eq.  (2-7)  is  over  polarization.  Eqs.  (2-4)-(2-7)  have  been 

derived  following  results  by  Wang  .  The  overlap  integrals  in  Eq.  (2- 
6)  account  for  the  mismatch  in  mode  profile  between  two  layers. 
Eqs.  (2-l)-(2-7)  are  used  repetitively  to  calculate  the  field  amplitude 
in  all  layers;  the  boundary  conditions  are  satisfied  only  for  a 
particular  set  of  values  of  gQ,X,  Bj,  and  An+2  (n  is  the  number  of 
layers). 

We  apply  this  model  to  a  particular  longitudinal  mode  (at  882.6 
nm)  of  the  SEL  structure  of  Fig.  2.  Fig.  15  shows  the  forward  and 
backward  output  power  versus  pumping  calculated  for  the  design 
longitudinal  mode  (at  882.6  nm)  of  the  structure  of  Fig.  2.  The 
quarter-wave  phase  shift  layer  effectively  diverts  most  of  the  power 
out  of  the  wafer.  It  is  seen  from  Fig.  16  that  the  spatial  intensity  is 
maximized  at  the  phase  shift  layer,  and  the  gain  is  saturated 

primarily  in  the  front  portion  of  the  resonator.  The  gain  in  the  back 
portion  of  the  resonator  is  available  for  extraction  by  the  next  mode 
with  shorter  wavelength  (at  860.8  nm),  since  its  intensity  is 

maximized  in  that  section  of  the  structure.  For  this  laser  structure 

the  threshold  gain  of  the  mode  at  860.8  nm  is  roughly  four  times 

that  of  the  main  mode  at  882.6  nm;  therefore  its  contribution  to  the 
output  spectrum  will  be  relatively  small  and  likely  be  as  a  broad 
amplified  spontaneous  emission  background.  In  fact  the 
experimental  results  obtained  from  this  structure  (cf.  Fig.  10)  seem 
to  agree  with  this  prediction. 
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SiTY  NORMALIZED  TO  SATURATION  INTENSITY 


0123456789  10  11 


Figure  16.  Spatial  intensity  distribution  for  (a)  the  Bragg  mode,  and 
(b)  next  higher  order  axial  mode 
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CHAPTER  4  SUMMARY  AND  DISCUSSION 


In  summary,  we  have  conducted  an  effort  in  the  development  and 
investigation  of  SEL  diodes  and  SEL  arrays,  and  have  made  major 
progress  and  accomplishments  during  this  18-month  project.  We 
developed  a  novel  SEL  using  a  transverse  buried  heterojunction 
structure  and  a  vertical  distributed  feedback  cavity,  and  successfully 
fabricated  and  demonstrated  the  operation  of  such  a  SEL.  The  SEL  has 
unique  characteristics  including  a  very  low  threshold  current  and  a  short 
optical  cavity  that  can  be  modulated  at  a  very  high  rate,  making  it 
important  for  optical  communications  and  signal  processing. 

We  further  extended  the  single  SEL  concept  to  the  formation  of  a  2-D 
array,  and  fabricated  and  demonstrated  a  SEL  array  with  6  x  16 
elements.  This  effort  proves  the  concept  that  this  SEL  design 
configuration  can  be  scaled  into  a  2-D  array  easily  and  has  a  potential  of 
forming  a  large  2-D  array  with  very  high  packing  density  on  a  wafer 
level.  These  features  are  very  desirable  for  optical  interconnect,  optical 
pumping,  and  signal  processing. 

We  also  undertook  a  theoretical  modeling  effort  in  conjunction  with 
the  experimental  development  effort,  with  the  objective  of  analyzing  and 
understanding  the  device  performance  characteristics  as  well  as 
providing  a  design  tool  and  guideline.  Two  theoretical  models  were 
developed.  One  model  is  based  on  the  Bloch  wave  approximation,  and 
another  model  uses  the  thin  film  discrete  wave  propagation  method. 
Both  models  have  proven  to  be  very  useful. 

During  the  course  of  this  project  several  technical  papers  were 
generated.  Four  papers  were  presented  at  the  technical  conferences  and 
two  papers  are  submitted  for  publication  in  refereed  journals.  The 
reprints  and  preprints  of  these  papers  are  attached  in  the  appendix. 

The  progress  and  achievements  we  have  achieved  during  this  project 
period  are  very  encouraging.  However,  these  results  are  preliminary, 
and  there  remains  a  number  of  technical  issues  regarding  the 
development  and  understanding  of  the  SEL  and  its  array  that  need  to  be 
addressed. 
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These  issues  include  (1)  the  relatively  low  quantum  efficiency,  (2)  the 
spatial  hole  burning  of  DFB,  (3)  the  broad  emission  linewidth,  (4)  the 
emission  beam  quality,  (5)  the  array  coherence,  (6)  the  power  scalability, 
and  (7)  the  development  of  DBR  type  SEL.  All  these  technical  issues  are 
very  important  in  understanding  the  device  physics  and  performance 
characteristics  of  the  SEL. 


CHAPTER  5  RECOMMENDATION 


Based  on  our  preliminary  findings  obtained  thus  far  during  the  course 
of  this  project,  we  found  the  SEL  technology  we  developed  has  a  number 
of  desirable  features  and  key  functional  advantages  over  the 
conventional  edge  emitting  laser  diodes  and  other  types  of  SELs. 
However,  because  of  the  remaining  technical  issues  we  mentioned 
previously,  the  ultimate  superiority  and  potential  limitations  of  this  SEL 
technology  cannot  be  candidly  assessed.  Therefore,  we  strongly 
recommend  that  a  follow-on  research  effort  in  this  SEL  technology  be 
continued.  The  follow-on  effort  will  investigate  those  technical  issues 
thoroughly  and  bring  this  technology  to  a  more  mature  stage  for  further 
assessment.  Such  an  effort  could  also  make  the  SEL  an  enabling 
technology  for  many  potential  DoD  applications  in  the  areas  of  optical 
signal  processing,  optical  communication,  optical  interconnect,  optical 
pumping,  and  optical  computing. 
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ABSTRACT 

A  threshold  current  of  2  mA  at  room  temperature 
CW  operation  is  realized  in  a  vertical  distributed  feedback 
surface  emitting  laser  diode  (VDFB-SELD)  with 
transverse  junction  buried  heterostructure  (TJBH).  In  this 
T J HI  1  structure,  the  vertical  distributed  feedback  active 
region  (AlGaAs/GaAs  multilayer)  is  entirely  surrounded 
by  N-  and  P-type  AlGaAs  cladding  layer  for  minority 
carrier  confinement.  The  far  field  angle  is  7  to  8  degrees. 
The  beam  shape  is  nearly  circular.  However,  the  lasing 
spectrum  is  broad  (  1.5  to  3  nm  )  compared  with  the  con¬ 
ventional  edge-emitting  laser.  Theoretical  model  of  the 
VDFB-SELD  TJBH  by  using  DFB  theory  and  its  computer 
simulation  of  the  spectrum  shows  good  agreement  with 
the  experimental  measurements. 


INTRODUCTION 

The  vertical  distributed  feedback  surface  emitting 
laser  diode  (VDFB-SELD)  is  very  advantageous  for  opto¬ 
electronic  integration  because  it  does  not  need  cleavage  nor 
a  backside  substrate  etching  process  to  form  a  vertical 
optical  cavity  [1,2].  High  reflectivity  of  the  cavity  mir¬ 
rors  and  high  efficiency  in  the  carrier  confinement  struc¬ 
ture  art  important  to  realize  low  threshold  current  opera¬ 
tion  of  the  VDI  B-SLLD.  High  reflectivity  of  the  cavity 
mirrors  was  already  realized  by  the  use  of  a  quarter 
wavelength  stack  of  AlGaAs  and  GaAs  thin  layers  [3]; 
however,  earner  confinement  was  not  adequate  in  the  pre¬ 
vious  structures  [1,2].  In  this  report,  we  realized  a  novel 
transverse  junction  buried  heterostnicture  (TJBH)  by 
selective  liquid  phase  epitaxy  (LPE)  and  selective  zinc 
diffusion,  and  achieved  a  room-  temperature,  low- 
threshold  current,  CW  operation  of  the  VDFB-SELD. 

DEVICE  FABRICATION 

Figure  I  depicts  the  VDFB  SELD-TJBH  structure. 
The  active  layer  is  made  of  a  quarter  wavelength  slack  of 
Ali,jGa,l7As  and  GaAs  thin  layers  prepared  by  metal- 
organic  chemical  vapor  deposition  (MOCVD).  The 
thicknesses  of  each  pair  of  AlGaAs/GaAs  layers  are  first 
chosen  and  then  adjusted  so  that  the  reflectivity  is  peaked 
near  the  gain  maximum  of  GaAs.  The  optical  cavity, 
formed  by  these  AlGaAs/GaAs  multilayers,  consists  of  20 
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pairs  at  the  top,  60  pairs  at  the  bottom,  and  a  phase 
shifler  GaAs  layer  in  between.  The  thickness  of  the  GaAs 
phase  shifting  layer  is  about  half  wavelength  of  the  gam 
maximum,  which  allows  the  laser  to  lase  at  the  Bragg 
wavelength  ol  the  multilayer.  This  multilayer  is  etched 
by  a  wet  etching  solution  leaving  mesas  with  different 
dimensions  ranging  from  2  x  8  to  3  x  17  pirn  on  the  top 
surface  and  10  j*m  in  height  with  the  DFB  type  active 
region.  The  N-typc  Al0uGa06As  cladding  layer  is  formed  I 

by  selective  LPE  using  silicon  nitride  on  top  of  the  mesa  ! 

as  a  selective  epitaxial  mask.  The  lateral  pn-junction  is 
formed  by  selective  zinc  diffusion  through  an  opening  in 
the  silicon  nitride  film  at  a  distance  of  3  to  4  ^m  away 
from  the  edge  of  the  mesa.  Due  to  lateral  diffusion,  the 
zinc  diffusion  front  is  inside  the  multilayer.  In  this  struc¬ 
ture,  the  GaAs/AlGaAs  part  of  the  multilayer  is  com¬ 
pletely  surrounded  by  the  N-  or  P-type  AlGaAs.  There¬ 
fore,  carrier  confinement  should  be  comparable  with  that 
of  buried  hetero-structures  used  wilh  the  edge  emitting 
laser  diode  .  The  zinc  diffusion  front  spreads  along  the 
cladding  region  as  well  as  along  the  AlGaAs/GaAs  mul¬ 
tilayer.  However,  carriers  are  injecled  predominantly  into 
the  the  GaAs  part  of  the  multilayer  because  the  turn-on 
voltage  for  an  AlGaAs  pn-junction  is  higher  than  for  a 
GaAs  pn  junction. 

DEVICE  RESULTS 

Figure  2  shows  the  light  output  versus  dc  current 
fl.-I)  characteristics  of  the  VDFB-SELD.  The  laser  is 
operated  CW  at  room  temperalure.  The  threshold  current 
is  2  mA  as  seen  in  the  inset  of  this  figure.  I  he  total  light 
output  increases  linearly  up  to  25  mA.  The  output  power 
is  28  jxW  at  10  mA  and  87  MW  at  25  mA. 

Figure  3(a)  shows  the  near-field  pattp"i  of  the 
VDFB-SELD.  The  light-emitting  region  is  tightly  confined 
within  the  rectangular  top  emission  surface.  Therefore, 
carrier  confinement  by  the  TJBH  structure  is  very 
effective.  Figure  3(b)  shows  the  far-field  pattern  at  an 
operating  current  of  20  mA.  It  shows  a  nearly  circular 
pattern  with  a  gaussian-like  intensity  profile.  The 
estimated  beam  divergence  angle  is  about  7  to  8  degrees 
This  indicates  a  good  spatial  coherence  in  this  device. 
However,  the  circular  emission  pattern  is  not  expected 
from  a  trapezoidal  mesa  cavity  with  a  rectangular  emis¬ 
sion  surface.  We  think  that  the  curved  boundary  between 
the  AlGaAs/GaAs  multilayer  active  region  and  AlGaAs 
cladding  region  may  affect  the  emission  beam  pattern. 

Figure  4  shows  the  CW  electroluminescence  specira 
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ol  the  VDI  B  SF.I  D- FJBII  structure. i  As  is  seen  in  the 
inset  ol  ll.i  figiuv.  I  he  tinsel  ol  sh.irperins.su>  n  lineal  884 
nn,  lx. nines  obvious  at  the  drive  elirrent  ol  I  .(>  iti A.  I  his 
tuireiil  level  is  consistent  woh  the  threshold  curient 
deduced  trom  the  1  1  characteristics.  The  hall  width  ol 
the  lasing  spectrum  at  884  nm  is  rather  wide  (  1.5  to  3 
nm  I  compared  with  an  edge  emitting  (I  I:)  laser  duxte 
anti  there  is  a  large  amplified  spontaneous  emission  at  the 
shorter  wavelength  side.  However,  there  are  fundamental 
differences  between  the  short -cavity  SI:  and  long-cavity 
1:1-1  Iasi  rs.  First.  the  mode  spacing  in  the  SFI  laser  is  much 
wider  than  the  1:1:  laser.  Tor  such  wide  sjiectral  sepaia- 
tion,  the  mode  briiadentng  mechanism  will  be  different 
from  that  ol  the  T.h  laser.  Second,  the  spontaneous  emis¬ 
sion  I  del  or  in  the  short  cavity  SI:  laser  is  expected  to  be 
much  larger  than  that  in  the  long  cavity  IT:  laser  [4). 
Both  these  (actors  will  affect  the  behavior  of  the  two 
types  of  lasers  in  terms  of  lasing  spectrum. 

t 

THEORETICAL  MODKI  ING  OK  DEVICE 
We  use  Bhxh  wave  theory  lor  the  DI  B  multilayer 
section  |5]  arid  consider  the  laser  as  an  amplifier  driven  by 
spontaneous  emission.  Figure  5  shows  the  idealized  struc¬ 
ture  we  used  for  the  calculation  ol  the  reflectivity  curve 
and  emission  spectrum.  E,  &  E4  ,  E2  &  E3  are  the  electric 
fields  of  the  forward-  and  backward  -going  waves  respec¬ 
tively.  at  the  two  interfaces  ol  the  phase  shiltcr  region 
( I  j ).  S,'s  are  the  equivalent  sources  at  the  interf  aces  due  to 
the  random  spontaneous  emission  in  each  section.  R.'s  and 
T,'s  are  the  reflection  and  transmission  coefficients  of  the 
Bloch  u-aves  at  the  interfaces.  R2.  R3  and  T2,  T3  refer  to 
the  phase  shifter  as  a  whole.  The  round-trip  condition  is 
given  by: 

(  1  — R,  R,  P/  )(  1-Rj  R4  P42  )  =  R,  R4  T2  Tj  P22  P42  (  I  ) 


where  the  P,  s  are  the  propagation  factors  in  each  section 
[5].  From  Tqn.f  1 ).  we  can  calculate  the  wavelengths  and 
threshold  gains  of  the  modes.  The  output  spectra!  power 
density  is  given  by: 
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where 

N  =  R2+  R4  P42  (  T2  T3  -  R2  R3  ) 

D  =  (  1  -  R,  R2  P/  )(  1  -  R3  R4  P4  )  -  R,  R4  T2  T3  P/  P42 

and  the  P.  (\)'s  are  the  spectral  power  densities  of  the 
equivalent  spontaneous  emission  sources  S,'s  at  the  four 
interfaces.  Realistic  gain  and  spontaneous  emission  spec¬ 
trum  were  used  for  the  calculations. 

COMPUTER  SIMULATION  RESULTS 

Figure  b  shows  the  reflectivity  curve  measured  from 
the  top  surface  and  the  theoretical  calculation  based  on 
Bloch  wave  theory.  The  agreement  is  very  good  in  predict¬ 
ing  the  shape  of  reflectivity  curve.  Figure  7  shows  the 
experimental  pulsed  spectrum  and  the  computer  calcu¬ 
lated  power  spectrum  from  F.qn.f  I )  and  l:qn.(2)  above,  at 


different  carrier  densities.  The  agreement  is  fairly  good. 

1  he  lowest  effective  threshold  power  gain  inode  ( C »,,,  —  i.t> 
cm  ')  Is  at  884  nrn  from  calculation,  all  the  other  modes 
w  hich  appear  as  the  humps  in  the  spectrum  have  a  much  | 
higher  effective  threshold  power  gain  (G,b—  l-2x  103cm  1  I. 

I  hus.  only  the  mode  at  884  nm  is  lasing  (i.e.  above  thres¬ 
hold).  all  the  other  modes  are  still  below  threshold  at 
these  current  levels.  Therefore,  the  emission  spectrum 
shows  a  sharp  peak  at  884  nm  and  much  broader  humps 
at  shorter  wavelengths  due  to  the  amplified  spontaneous 
emission  prexiuced  by  these  high  threshold  gam  modes.  In 
this  model,  we  assume  that  the  earner  density  is  given 
and  the  gain  spectrum  is  then  calculated  with  a  semi- 
empirical  model,  assuming  Fermi -Dirac  distributions  in 
the  bands.  Gain  saturation  is  not  taken  into  account. 

1  herefore.  the  model  is  really  valid  only  below  threshold 
and  fails  to  give  the  exact  amplitude  of  the  lasing  peak 
above  threshold.  The  fact  that  it  reproduces  the  shape  of 
the  spectrum  above  threshold  elsewhere  indicates  that 
there  is  inhomogeneous  broadening  and  that  the  carrier 
density  is  not  damped  at  threshold.  The  measured 
I  WHM  of  the  lasing  peak  ts  probably  larger  than  the  real 
one  because  of  the  limited  resolution  (7  A )  of  the  meas¬ 
urement  (  which  was  incorporated  in  the  simulations). 

CONCLUSIONS 

In  conclusion,  the  VDFB-SELD-TJBH  structure  is 
realised  with  the  combination  of  MOCVD.  selective  l.PL, 
and  zinc  diffusion  techniques.  In  this  structure  ,  an 
MOCVD-grown  AlGaAs/GaAs  multilayer  is  surrounded 
completely  by  P-  and  N-type  AlGaAs  cladding  layers, 
such  that  lateral  earner  injection  and  carrier  confinement 
are  achieved.  Nominal  threshold  current  for  CW  operation 
Is  2  mA.  The  radiation  region  is  well  confined  within  the 
active  region.  The  far-held  pattern  is  circular  and  well  I 
defined.  However,  lasing  spectra  at  884  nm  is  broad  ( 1.5 
to  3  nm)  and  ihere  is  large  amplified  spontaneous  emis¬ 
sion  at  shorier  wavelengths.  Compuler  simulation  shows 
good  agreement  with  the  measurement.  Furl  her  improve¬ 
ment  of  the  spectrum  can  be  realized  by  propter 
modifications  of  the  opiical  cavity  and  layer  thickness  of 
the  multilayer  and  phase  shifter. 
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FIGURE  CAPTIONS 

(a  (Schematic  diagram  of  the  VDFB-SELD-TJRH 
structure.  ( b)Cross-sectional  SEM  picture  of  the 
VDEB-SELD-TJBII  structure.  The  white  con¬ 
trasted  region  at  the  right  side  of  the  multilayer 
is  a  zinc  diff used  P-type  AlGaAs  cladding  layer. 
Light  output  versus  dc  drive  current  of  the 
VDFB-SELD-TJBH  structure. 

(a)Near-field  pattern  of  the  VDFB-SEI.D.  Also 
shows  the  line  profile  of  the  field  intensity 
across  the  long  edge  of  the  rectangular  emission 
region.  (b)Far-field  pattern  and  its  intensity 
line  profile  across  the  beam  center  of  the  VDFB- 
SE1.D. 

CW  electroluminescence  spectra  of  the  VDFB- 
SEI.D. 

Schematic  diagram  of  the  Idealized  VDFB-SELD 
structure  used  for  theoretical  model.  L3  is  the 
phase  shifter  region.  L2&L4&LS  are  the  mul¬ 
tilayer  DFB  section. 

Measured  and  calculated  reflectivity  curve  of  the 
multilayer  from  the  top  surface. 

Spectrum  of  the  VDFB-SELD.  (a (Experimental 
measurements  (400  ns  pulse).  Resolution  of  the 
monochrometer  is  0.7  nrn  in  the  measurement, 
(b (Computer  simulation  results  for  carrier  den¬ 
sity  N=-0.5xlO'*  cm  3  ( I ),  IxlO1*  cm  3  (2).  2xl018 
cm  1  (3).  3xlOls  cm  3  (4). 
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tition  rate  can  only  be  achieved  In  highly  doped 
(50%)  crystals  where  by  4fn  2  *  4/,s  z.  4/«3  7  * 

%  2  cross  relaxation  the  long- living  lower  laser 
level  can  be  depopulated  Highly  doped  flash- 
lamp-pumped  YAIO3  Er  has  a  typical  laser  thresh¬ 
old  of  -30  J  for  a  crystal  of  75-mm  length  and 
6  35-mm  diameter  100-Hz  operation  would  thus 
lead  to  an  average  pump  power  exceeding  3  kW. 
well  beyond  the  thermal  damage  threshold 
Therefore,  a  reduction  of  the  threshold  is  required 

As  most  of  the  pump  light  of  a  xenon  flashlamp 
is  absorbed  far  above  the  upper  laser  level,  this 
excitation  leads  to  cross  relaxations,  emission  of 
resonance  lines,  and  stimulated  muftiphonon  tran¬ 
sitions  These  processes  can  enhance  the  laser 
threshold  or  cause  crystal  heating.  To  reduce  the 
threshold,  the  pump  light  was  spectrally  filtered  by 
introducing  color  filters  or  dyes  In  the  pump  cavity 
between  the  flashlamp  and  the  laser  rod  Experi¬ 
mentally  a  reduction  of  the  absorbed  energy  In  the 
crystal  at  laser  threshold  by  a  factor  of  ^4  and  a 
pulse  repetition  rate  up  to  125  Hz  could  be 
achieved  (12  min) 

MJ3  Optimization  of  th«  gain  In  Nd:SHica  single- 
mode  fibers 

M  DIGONNET.  K  LIU.  H.  J  SHAW.  Stanford  U., 
Ginzton  Laboratory.  Stanford.  CA  94305-4085. 

The  development  of  efficient  (low  pump  power) 
optical  amplifiers  and  sources  based  on  rare-earth 
doped  or  other  active  fibers  depends  critically  on 
the  optical  gain  of  the  fiber,  especially  when  three- 
level  or  low  gain  (e  g  ,  nonlinear)  processes  are 
involved  In  optically  end-pumped  fibers,  the  gain 
depends  on  optical  energy  confinement  provided 
by  the  fiber  and  consequently  on  the  fiber  N.A  and 
core  dimension  These  dependences  were  Inves¬ 
tigated  experimentally  and  theoretically  to  estab¬ 
lish  the  range  of  fiber  parameters  which  optimizes 
the  gain.  In  Nd:SiO?  fibers  from  different  manu¬ 
facturers.  we  have  measured  an  optical  gain  rang¬ 
ing  from  0.15  to  0  4  dB/mW.  These  data  were 
analyzed  with  a  model  of  gain  in  fibers  accounting 
for  energy  confinement  through  pump  and  signal 
mode  spatial  overlap  It  shows  that  the  observed 
gain  differences  are  largely  due  to  differences  In 
the  fiber  Index  profiles  as  opposed  to  differences 
in  doping  levels.  This  suggests  that  the  fiber  gain 
can  be  optimized  by  selecting  the  core  size  and 
N.A  in  a  well-defined  range,  predicted  theoretical¬ 
ly  and  supported  by  gain  measurements.  As  much 
as  a  2  -  4  gain  increase  (and  equivalent  reduction  In 
the  required  pump  power)  over  what  we  have  ob¬ 
served  may  be  expected  from  careful  fiber  design. 

(12  min) 

MJ4  Effect  of  ion- Ion  Interactions  on  the  life¬ 
time*  of  the  29 00 -nm  laser  levels  of  Er3* 

W.  Q  SHI.  J.  MACHAN.  M.  BIRNBAUM.  U.  South¬ 
ern  California,  Center  for  Laser  Studies,  Los  Ange¬ 
les.  CA  90089- 1 1 12;  M.  BASS.  U.  Central  Florida, 
Orlando  PL  327 16-000 1 ;  M  KOKTA,  Union  Carbide 
Corp  .  Washougal.  WA  98671. 

The  effects  of  Ion-ion  interactions  In  heavily 
doped  ErYAG.  Er,Nd:YAG,  and  Er.Nd.YALO  on  the 
lifetimes  of  the  4/t1,2  and  Vn,2  levels  of  Er3*  were 
studied.  In  Er  YAG,  the  4/13,?  lifetime  first  In¬ 
creases  with  increasing  concentration  and  then 
decreases  The  Increase  can  be  attributed  to  radi¬ 
ation  trapping  and  the  decrease  to  concentration 
quenching  In  Nd(1%),  Erf  15%)  YAG.  the  4/,„2 
level  lifetime  is  halved,  while  the  4(13 2  ievel  life¬ 
time  suffers  an  approximately  sixty-fold  reduction, 
but  Tn,2  >  fn/?-  However,  study  of  a  YALO 
crystal  (1%  Nd.  15%  Er)  disclosed  a  r,,  ?/rI3,2 
ratio  of  ~  1.7.  Indicating  the  feasibility  of  four. level 
laser  operation  Ion-Ion  Interaction  coefficients 
were  obtained  from  an  analysis  using  the  rata 


equation  model  Laser  operation  at  2900  nm  with 
Nd.Er  YAG.  and  Nd.ErY  ALO  is  presented 

(12  min) 

MJ5  Semiconductor  laser  theory  with  many- 
body  effects 

STEPHAN  W  KOCH.  HARTMUT  HAUG,  MURRAY 
SARGENT  III.  U  Arizona.  Optical  Sciences  Center. 
Tucson.  A Z  85721 

A  theory  of  a  semiconductor  laser  is  developed 
that  includes  the  many-body  effects  due  to  Cou¬ 
lomb  interactions  The  theory  Is  valid  tor  both  3-D 
bulk  semiconductors  as  well  as  quasi- 2-D  quantum 
well  structures  We  emphasize  plasma  density- 
dependent  band  gap  renormalization,  broadening 
due  to  intraband  scattering,  and  electron  hole 
Coulomb  enhancement  The  very  short  intraband 
scattering  relaxation  time  allows  us  to  eliminate 
the  interband  polarization  adiabatically  and  to  in¬ 
troduce  a  hydrodynamic  description  of  the  intra¬ 
band  kinetics  From  this  general  formulation  a 
diffusion  equation  for  the  carrier  density  is  derived. 
The  resulting  diffusion  coefficient  decreases  with 
carrier  density  and  laser  intensity  due  to  the  reduc¬ 
tion  of  the  electron  drift.  We  use  our  theory  in  the 
problem  of  laser  gain,  index,  and  side  mode  insta¬ 
bilities  We  show  that  near  the  laser  operating 
point,  our  many-body  theory  can  be  approximated 
by  a  simple  rate  equation  formalism  However,  in 
contrast  to  the  usual  rate-equation  theory,  the 
semiconductor  rate  constants  are  functions  of 
temperature,  tuning,  and  carrier  density  ( 12  min) 

MJ6  Noise  characteristics  of  nearly  single¬ 
mode  semiconductor  lasers 

G.  GRAY.  RAJARSHI  ROY,  Georgia  Institute  of 
Technology.  School  of  Physics.  Atlanta,  GA 
30332. 

We  have  developed  a  computer  model  which 
simulates  the  nonlinear  Langevin  equations  for  a 
semiconductor  laser  with  one  strong  mode  and  one 
small  side  mode  With  this  model  we  are  able  to 
probe  the  laser  dynamics  over  a  wide  range  of 
operating  powers,  from  below  threshold,  when  the 
laser  Is  essentially  a  thermal  noise  generator,  to 
way  above  threshold.  The  model  correctly  repro¬ 
duces  the  well-known  intensity  noise  spectra  and 
mode  partition  effects.  We  further  calculate  noise 
variances  and  intensity  autocorrelation  and  cross¬ 
correlation  functions  for  both  modes  as  weft  as  for 
the  total  power,  and  we  compare  these  results  to 
linearized  theories.  Finally,  we  calculate  the  laser 
field  spectrum  and  investigate  the  effect  of  the  side 
mode  on  the  laser  Mnewidth.  All  the  above  proper¬ 
ties  are  calculated  for  many  different  values  of  the 
mode  suppression  ratio.  (12  min) 

MJ7  Kinetics  of  gradual  degradation  In  semi¬ 
conductor  lasers 

YU  L.  KHAIT  AL,  Technion  -  Israel  Institute  of  Tech¬ 
nology.  Haifa  32000.  Israel;  JOSEPH  SALZMAN, 
Technion -Israel  Institute  of  Technology.  Solid 
State  Institute.  Haifa  32000.  Israel;  R  BESERMAN, 
Technion-  Israel  Institute  of  Technology,  Depart¬ 
ment  of  Electrical  Engineering.  Haifa.  32000.  Isra¬ 
el. 


bE  is  localized  around  this  a'om  In  the  present 
model,  such  synchronization  is  inherent  the  fluc¬ 
tuating  atom  forms  a  transient  point  defect  ebich 
behaves  as  a  nonradiative  rero'iihmaiion  center  r 
lifetime  At  ^10  "  10  s  f  1  e*»  c ar r ier s  al 

trapped  at  these  centers  leading  to  synchronoA 
release  of  localized  energy  E  this  model  is 
based  on  die  kinetic  many  body  theory  of  s horf- 
tived  large  energy  fluctuations  of  atoms  *  The 
laser  degradation  rate  K  is  given  by  K  -  Kc 
exp|  -(E,  -  bE)/kT\  with  K  and  <>£.  expressed  in 
terms  ot  the  laser  parame'eis  ttarrier  concentra¬ 
tion.  energy  (jap.  detect  oonceitti.ition.  etc  )  Cal¬ 
culated  laser  lifetimes  are  in  good  agreement  with 
experimental  data  ( 12  min) 


1  J  D  Weeks.  J  C  Fully,  and  L  C  Ktmerlmg. 
Phys  Rev  B  12.  3286  ( 19-’5\ 

2  Yu  L  Khait.  Phys  Rep  99.  237(1983) 

3  Yu  L  Khait  and  R  Be sor man  Phys  Rev  B  33, 
2983  (1986) 


MJ8  Degradation  of  mode  discrimination  in  dis¬ 
tributed  feedback  lasers  due  to  spatial  hole  burn¬ 
ing 

WILLIAM  S  RABINOVICH.  B  J  FELDMAN.  tj  S 
Naval  Research  Laboratory.  Laser  Physics  Branch, 
Washington.  DC  20375 

Distributed  feedback  (DEB)  lasers  have  become 
of  Increasing  interest  in  fields  which  require  single 
longitudinal  mode  (SLM)  sources.  This  SLM  be¬ 
havior  occurs  because  the  threshold  gam  for  a 
longitudinal  mode  in  a  DFB  laser  increases  as  that 
mode's  frequency  moves  away  from  the  Bragg 
frequency.  We  have  derived  coupled-mode  equa¬ 
tions  for  the  above  threshold  behavior  of  homoge¬ 
neously  broadened  DFB  lasers,  including  the  spa¬ 
tial  hole  burning  grating  induced  by  the  lasefl 
modes.  These  equations  were  numerically  inteH 
grated  for  boundary  conditions  appropriate  to  a 
variety  of  DFB  structures  To  study  the  mode 
discrimination  of  the  DFB  laser,  we  allow  one  lon¬ 
gitudinal  mode  to  exceed  threshold  and  induce  a 
spatial  hole  burning  grating  We  simultaneously 
solve  for  the  threshold  of  a  second  longitudinal 
mode  in  the  presence  of  both  this  field-induced 
grating  and  the  externally  Induced  grating  The 
ratio  of  the  threshold  of  this  second  longitudinal 
mode  to  the  average  saturated  line  center  gain 
available  In  the  laser  is  defined  as  the  mode  dis¬ 
crimination.  The  higher  the  discrimination  the 
more  likely  the  laser  is  to  be  SLM  Our  analysis 
shows  that  the  mode  discrimination  of  DFB  lasers 
with  nonsaturable  Index  coupling,  such  as  semi¬ 
conductor  DFB  lasers,  is  seriously  degraded  as  the 
*  laser  exceeds  output  powers  a  few  times  above 
threshold.  The  same  does  not  hold  true  for  gain 
coupled  DFB  lasers.  Indicating  that  these  types  of 
laser  may  be  more  suitable  for  high  output  power 
SLM  applications.  ( 12  min) 


MJ9  Single-mode  model  of  surface-emitting 
DFB  and  DBR  semiconductor  lasers 


THADOEUS  G.  D2IURA.  S  C  WANG.  Lockheed 
Palo  Alto  Research  Laboratory,  Palo  Alto.  CA 
94304 


Gradual  degradation  in  semiconductor  lasers  is 
caused  by  the  creation  and  migration  of  point  de¬ 
fects  which  require  jumping  of  atoms  from  their 
initial  position  to  disordered  ones  (with  activation 
energy  E,  of  a  few  electron  volts)  The  injection  of 
a  nonequilibrium  electron  hole  plasma  may  in¬ 
crease  the  probability  of  structural  changes  and 
reduce  their  activation  energy  by  bE  ~~  Fr  —  F„ 
through  nonradiative  recombination'  (F  and  Fx  are 
the  quasi-Fermi  levels)  Such  enhancement  Is 
possible  only  If  the  recombination  is  synchronized 
with  the  atom  jump  event  and  the  released  energy 


The  surface-emitting  semiconductor  laser1  3  is 
expected  to  he  a  key  component  in  future  opto¬ 
electronic  and  optical  signal  processing  applica¬ 
tions  Several  designs  a»®  under  development 
Including  electrically  excited  and  phoiopumped  lag 
sers  The  requirement  of  a  surface-emitting  qe^ 
ometry  imposes  important  constraints  on  laser  de¬ 
sign.  which  must  be  considered  to  achieve  an 
efficient  emitter  For  e»anip)*».  the  method  of 
extracting  the  laser  heam  must  acr  ommndate  that 
part  of  the  structure  providing  *»tecti  k  il  *'»citat'on 
Also,  a  model  that  calculates  laser  proppHies  only 
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at  threshold  is  clearly  not  adequate  We  have 
developed  a  transfer  equation  model  of  short  cavi¬ 
ty  OFB  and  DBR  surface-emitting  Users  in  which 
the  resonator  is  defined  by  multilayer  quarter- 
wavelength  semiconductor  structures  providing 
distributed  feedback  These  lasers  naturally  lend 
themseUes  to  a  description  in  which  the  forward 
and  backward  waves  are  propagated  from  layer  to 
layer  by  a  set  of  equations  detailing  the  gam  and 
loss  due  to  the  material  and  the  feedback  and 
radiation  loss  due  to  the  Bragg  grating  We  apply 
this  method  to  calculate  power  output,  longitudinal 
modes  and  axial  intensity  profiles  for  several  laser 
types  In  particular  we  show  that  the  forward  to 
backward  power  ratio  can  depend  on  mode  num¬ 
ber  as  well  as  pumping  level  ( 12  min) 

1.  M  Ogura.  W  Hsin.  M  C  Wu,  S  Wang,  J.  R 
Whmnery.  and  S  C  Wang.  Appi  Phys  Lett  51, 
1666  (1987) 

2  F  Koyama.  ft  Tomomatsu.  and  K  Iga.  in  Tech¬ 
nical  Digest.  Conference  on  Lasers  and  Electro- 
Optics  (Optical  Society  of  America.  Washing¬ 
ton,  DC.  1988),  paper  W05. 

3  Z  L  Uau,  J  N  Walpole,  and  V  Diadiuk.  in 
Technical  Digest.  Topical  Meeting  on  Integrated 
and  Guided  Wave  Optics  (Optical  Society  of 
America.  Washington,  DC.  1988).  paper  WA1 
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MK  1  Detection  of  phase-modulated  optical  sig¬ 
nals  through  two-wave  mixing  In  photorefractlve 
materials 

L.  BOUTSIKARIS,  FRED  DAVIDSON.  Johns  Hop¬ 
kins  U .  Electrical  Engineering  Department.  Balti¬ 
more.  MD  21218 

Two- wave  mixing  in  photorefractlve  crystals 
has  been  shown  to  combine  coherently  a  modulat¬ 
ed  optical  signal  beam  with  a  strong  pump  optical 
field  1  For  a  phase-modulated  signal  beam,  the 
induced  index  of  refraction  phase  grating  inside  the 
crystal  yields  a  coupling  of  the  optical  fields  of  the 
pump  and  signal  so  that  subsequent  Intensity  de¬ 
tection  results  in  apparent  Intensity  modulation  of 
the  transmitted  signal  or  pump  beams.  This  pro¬ 
vides  an  easy  way  to  detect  phase  differences 
between  two  optical  beams.  Some  of  the  advan¬ 
tages  the  above  scheme  offers  over  commonly 
used  coherent  optical  phase  detection  schemes 
are;  no  loss  of  3  dB  in  the  signal  beam  intensity  at 
the  50/50  beam  splitter  (used  to  combine  signal 
and  local  oscillator  fields);  automatic  alignment  of 
the  index  of  refraction  grating  with  the  two  incident 
optical  wavefronts;  insensitivity  to  slow  random 
changes  in  the  phases  of  both  beams.  Experi¬ 
mental  results  are  given,  and  a  new  definition  for 
the  resulting  mixing  gain  for  the  phase-modulated 
beam  is  introduced,  through  which  optimal  condi¬ 
tions  for  input  beam  intensity  ratios  can  be  found 
Applications  to  coherent  homodyne  detection  are 
discussed  (12  min) 

1  F.  Davidson,  L.  Boutslkarls,  and  M.  Krainak, 
Coherent  Optical  Detection  Through  Two- 

Wave  Mixing  in  Photorefractlve  Materials, "  Opt 

Lett  13,500  (1988) 


MK2  Broadband  communication  system  using  a 
standing  wave  acoustooptic  modulator 

» 

R  S  TIAN.  CHARLES  S.  IH.U  Delaware,  Electrical 
Engineering  Department,  Newark,  DE  197  16. 

There  is  great  need  for  an  all-service  communi¬ 
cation  system  by  which  broadband  analog  and  high 
speed  digital  information  can  be  transmitted  simul¬ 
taneously  and  freely  Such  a  system  can  best  be 
implemented  with  frequency-division-multiplexmg 
(FDM).  The  effective  bandwidth  of  such  a  system 
is  in  the  multigigahertz  range  Because  of  the 
enormous  bandwidth,  the  optical  fiber  communica¬ 
tion  (OFC)  system  is  a  good  candidate  for  such 
applications  However,  because  of  the  serious 
intermodulation  problems,  it  is  still  difficult  to  pro¬ 
vide  such  services  using  OFC.  A  novel  method, 
now  known  as  double-beam-modulation  (DBM) 
was  suggested  and  demonstrated  using  a  standing- 
wave  acoustooptic  modulator  (SWAOM)  to  216 
MHz  The  DBM  can  be  implemented  in  many 
forms  '  We  describe  the  experiments  of  DBM 
using  SWAOM  for  frequencies  beyond  1  GHz. 
Theoretical  analyses  and  comparisons  of  DBM 
with  conventional  FDM  systems  are  described. 

(12  min) 

1.  C  S.  Ih.  R  G  Hunsperger.  J  J  Kramer,  R  Tlan, 
X.  Wang.  K.  Kissa.  and  J  Butler,  "A  Novel 
Modulation  System  for  Optical  Communica¬ 
tions."  Proc  Soc  Photo-Opt  Instrum  Eng.  878 
(1988). 

MK3  Novel  Integrated-optics  modulator  for  op¬ 
tical  communications 

KARL  KISSA.  X.  L.  WANG.  CHARLES  S.  IH,  R.  G. 
HUNSPERGER.  U.  Delaware,  Electrical  Engineer¬ 
ing  Department.  Newark.  DE  197  16. 

An  integrated-optics  standing-wave  surface- 
acoustic-wave  modulator  (SWSAWM)  for  optical 
communications  is  being  investigated.  This  novel 
modulation  method,  which  has  been  described  pre¬ 
viously,  ’  was  used  to  upconvert  a  color  TV  video 
signal  to  TV  channel  4  by  generating  a  70-MHz 
subcarrier  with  a  Ti-LiNbOj  modulator.  Prisms 
were  used  to  couple  light  in  and  out  of  the  wave¬ 
guide.  A  photodetector  placed  at  one  of  the  dif¬ 
fracted  beams  received  the  upconverted  signal. 
An  ordinary  TV  receiver  recovered  the  original 
video  signal  Picture  quality  was  good.  The 
phase  of  the  acoustooptlcally  generated  subcarrier 
could  be  changed  by  changing  the  phase  relation¬ 
ship  between  the  rf  signals  driving  the  two  inter- 
digital  transducers.  The  standing  wave  nature  of 
the  modulator  was  also  demonstraled  by  the  fre¬ 
quency  dependence  of  the  diffracted  optical  pow¬ 
er.  Nonlinearities  In  surface  acoustic  wave  gener¬ 
ation  and/or  acoustooptic  diffraction  caused  sub¬ 
carrier  distortion  on  high  diffraction  efficiencies. 
This  effect  is  being  examined  theoretically  and 
experimentally  to  determine  its  Influence  on  sys¬ 
tem  performance.  Subcarrier  generation  at  high¬ 
er  frequencies  (600  MHz  and  above)  is  also  being 
explored  (12  min) 

1.  C.  S.  Ih.  R  G.  Hunsperger,  J.  J.  Kramer,  R  Tlan, 
X  Wang.  K  Kissa.  and  J.  Butler,  A  Novel 
Modulation  System  for  Optical  Communica¬ 
tions."  Proc  Soc  Photo-Opt.  Instrum  Eng  876 
(1988) 

MK4  Use  of  edge  light  emitting  diodes  In  signal 
coherence  multiplexing 

THEODORE  L  HOUK,  R.  AARON  FALK.  RAYMOND 
W  HUGGINS.  Boeing  Co..  P  O  Box  3999,  Seattle. 
WA  98124-2499 

An  edge  light  emitting  diode  (ELED)  emission 


spectrum  is  near  Gaussian  at  820-nm  center  wave¬ 
length  with  33  nm  FWHM.  promising  a  short  coher¬ 
ence  length  within  the  range  of  feasible  titanium 
tndiffused  lithium  niobate  (Ti.LiNbOj  integrated  op¬ 
tics  unequal  arm  Mach-Zehnder  inter ferometer 
(IOMZ)  optical  path  differences  Interferometric 
ELED  coherence  measurements  were  made  in  a 
tilted  mirror  Michelson  interferometer.  Fringes 
were  viewed  with  a  near-lR  video  camera  and 
digifaf  imaging  frame  store  connected  to  a  comput¬ 
er  This  apparatus  measured  ail  the  quantities 
needed  to  compute  the  coherence  function  (Born 
and  Wolf),  which  fits  a  Gaussian  lunction  of  6  4-^m 
standard  deviation  (SO)  and  peak  value  0  78  We 
observed  interference  fringes  in  the  output  light  of 
the  second  of  two  Michelson  interferometers  in 
series  whenever  any  two  of  the  four  possible  opti¬ 
cal  path  lengths  were  equal  within  3  SDs.  20  pm 
These  fringes  are  the  reconstructed  signal  for  co¬ 
herence  multiplexing.  Proton  exchange  applied  to 
one  side  of  a  symmetric  ThLINbOa  IOMZ  can  form 
a  2000-pm  optical  path  difference  via  changes  m 
the  substrate  refractive  index  From  such  MZs. 
we  propose  to  form  coherence  multiplexing  sys¬ 
tems  of  up  to  ten  parallel  channels  ELED  power 
attenuated  by  system  T-junctions,  background,  and 
detector  noise  imposes  this  channel  limit. 

( 12  min) 

MK5  All-optical  demultiplexing  by  frequency 
upconversion 

HONG  LIN.  JOHN  STEPHEN  SMITH.  UC-Berkeley. 
Electrical  Engineering  &  Computer  Sciences. 
Berkeley.  CA  94720. 

We  propose  the  use  of  frequency  upconversion 
for  time  division  demultiplexing.  Suppose  we  re¬ 
ceive  an  N-channel  signal  at  frequency  f\  and 
make  a  local  short-pulse  laser  at  frequency  f 2 
have  the  same  repetition  rate  as  one  of  the  chan¬ 
nels.  Mixing  them  properly  In  a  nonlinear  crystal, 
we  can  upconvert  only  one  channel  to  frequency 
73  =*  M  +■  72.  We  can  demultiplex  N  channels  by 
N  such  units  The  frequency  upconversion  is  a 
very  fast  (femtosecond  domain  or  10E5  GHz)  and  a 
low  noise  process.  We  show  that  nearly  80% 
conversion  efficiency  of  1.3-0  58  ^m  can  be  ex¬ 
pected  for  50- W  peak  power  pumping  at  1 .06  pm 
in  a  3-mm  long  barium  sodium  niobate  crystal  1 
The  use  of  waveguide  structure  or  high  nonlinear 
organic  materials  may  relax  the  pumping  power 
requirement  substantially 

We  discuss  several  techniques  for  increasing 
efficiencies  of  the  system  and  the  possibility  of 
using  semiconductor  lasers  as  pumps  The  prob¬ 
lems  of  on-off  ratio  and  crosstalk  shoufd  not  exist 
with  a  properly  shaped  pulse  In  summary,  we 
show  that  the  all-optical  demultipler  by  frequency 
upconversion  Is  a  very  Interesting  device  for  a  high 
bit  rate  time-division-multiplexing  optical  commu¬ 
nication  system  2  An  experiment  is  under  way  in 
our  laboratory.  ( 12  mtn) 

1.  F.  Zernfke  and  J.  £  Midwinter,  Applied  Nonlin¬ 
ear  Optics  (Wiley,  New  York,  1973) 

2.  R  S  Tucker  et  a/  ,  in  Photonic  Switching,  T  K 

Gustafson  and  P  W  Smith.  Eds  (Springer-Ver- 

lag.  Berlin,  1988).  p.  206 

MK6  Free  apace  optical  pulse  position  modula¬ 
tion  communication  with  laser  diode  transmitter 
and  avalanche  photodiode  receiver 

X  SUN,  FREDERIC  OAVIOSON,  Johns  Hopkins  U  , 
Electrical  Engineering  Department.  Baltimore.  MD 
21218 

A  direct  detection  optical  communication  sys¬ 
tem  for  use  in  an  intersatellite  link  was  constructed 
with  an  AlGaAs  laser  diode  (Hitach-HL83  ME  \  = 
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Figure  1  Schematic  diagram  ot  the  transverse  buned  heteroiunction 
surface  etmuing  laser  ifcnir 
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ABSTRACT 

We  extend  the  Bloch  wave  calculation  method  to  the  case  of  a 
periodic  structure  where  the  gain  Goss)  and  the  index  of  refraction  depend 
.  on  frequency.  We  establish  a  procedure  to  find  the  propagation  constants 
and  the  reflection  and  transmission  coefficients  for  the  Bloch  waves  in  the 
general  case,  which  has  not  been  treated  previously.  We  then  use  the 
results  of  this  analysis  to  construct  a  model  for  a  surface-emitting  laser 
diode  (SELD)[2,3],  which  has  a  vertical  distributed  feedback  structure. 
This  model  considers  the  laser  as  an  amplifier  driven  by  spontaneous 
emission  to  compute  the  emission  spectrum.  At  the  present  time  it  is  valid 
only  below  threshold.  We  present  experimental  results  and  compare  the 
measured  emission  spectrum  with  the  theoretical  one.  The  reasons  for  the 
discrepancies  and  possible  improvements  to  the  device  and  to  the  theory 
are  then  discussed. 
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I.  Introduction 

Propagation  of  waves  in  periodic  structures  has  been  studied  extensively  in  the  past 
and  a  great  number  of  methods  has  been  used  to  solve  the  problem.  (For  a  good  sum¬ 
mary  and  references  on  the  subject,  see  the  review  paper  by  Elachi  [1].)  There  are  two 
main  types  of  such  structures  that  are  of  practical  interest  for  semiconductor  laser  diodes 
(Fig.  1):  (a)  periodic  waveguides  and  (b)  multilayers.  The  second  case  is  encountered  in 
surface-emitting  laser-diodes  (SELD)[2,3j.  These  are  basically  distributed  feedback 
(DFB)  semiconductor  lasers,  with  a  vertical  cavity  and  lateral  p-n  junction.  The  first 
order  grating  is  formed  by  alternating  layers  of  GaAs  and  GaAlAs.  This  structure  is  very 
advantageous  for  optoelectronic  integration,  since  it  does  not  need  cleavage  nor  backside 
etching  to  form  a  vertical  cavity,  as  some  other  S ELD’s  require  [4], 

In  most  of  the  previous  studies,  the  material  characteristics  of  the  periodic  structures 
were  supposed  to  be  constant  across  the  band  of  frequency  of  interest  This  is  a  good 
approximation  in  usual  edge-emitting  distributed  feedback  (DFB)  or  distributed  Bragg 

*  Electrotechnical  Laboratory  (L  I  L),  1-1-4  Umezono,  Sakura-mura  Niihah-gun,  Ibaraki  305,  Japan. 

**  Lockheed  Mitsilei  A  Space  Company,  Palo  Alto  Research  Lab,  325 1  Hanover  Street,  Palo  Alto,  CA  94304. 


reflector  (DBR)  laser- diodes,  where  we  are  looking  at  bands  that  are  2  to  3  nm  wide. 
However,  for  SELD’s,  because  of  the  strong  coupling  coefficients  required,  the  stop-band 
has  a  width  of  35  nm  or  more  'in  the  wavelength  region  of  850  to  90C  am).  When  we 
have  such  a  large  bandwidth,  material  dispersion  cannot  be  neglected  anymore.  The 
index  of  refraction  can  change  apprecianly,  but  the  main  effect  is  gain  and  loss  variations 
•with  frequency  in  the  semiconductor  material.  Also,  because  we  have  a  multilayer  where 
only  the  GaAs  is  pumped,  we  have  a  simultaneous  index  and  gain  grating. 

In  this  paper,  we  first  give  a  summary  of  the  modified  eigenmode  theory  that  takes 
these  effects  into  account.  We  also  consider  briefly  what  happens  when  the  periodic 
structure  does  not  have  an  even  sjTnmetry  anymore.  (Most  of  the  previous  studies  con¬ 
sidered  structures  with  even  symmetry  only.)  Then,  we  check  the  results  of  eigenmode 
theory  by  solving  exactly  the  two  layers  case.  Next,  we  derive  the  reflection  and 
transmission  coefficients  for  the  interface  between  a  periodic  waveguide  and  a  uniform 
waveguide.  Using  these  results,  we  are  able  to  compute  the  total  reflectivity  of  the  DFB 
structure  used  to  fabricate  the  SELD,  as  a  function  of  wavelength  and  we  compare  this 
result  with  the  measured  reflectivity.  We  then  summarize  the  experimental  data  of  the 
fabricated  device.  Next,  we  compute  the  theoretical  emission  spectrum  (below  thres¬ 
hold)  and  compare  it  to  the  measured  emission  spectrum.  The  model  used  for  this  calcu¬ 
lation  considers  the  laser  as  an  amplifier  driven  by  spontaneous  emission.  Finally,  we 
discuss  possible  improvements  to  the  theory  and  to  the  device. 

We  will  use  an  extension  of  the  eigenmode  theory  for  periodic  structures,  as  first 
developed  in  [5].  A  good  exposition  can  be  found  in  [6].  The  advantage  over  coupled 
mode  theory  (see  for  example  [?])  is  that  we  are  really  working  with  independent  modes, 
which  makes,  in  our  opinion,  the  boundary  conditions  much  easier  to  express.  (Although 
this  has  been  much  disputed,  see  [8,9].)  Another  advantage  is  that  it  allows  us  to  get 
analytical  formulas  for  the  output  spectrum  of  a  DFB  laser.  This  is  not  possible  with 
coupled  mode  theory  [10]. 
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IL  General  Theory 

The  equation  for  the  propagation  of  a  scalar  wave  in  a  one-dimensional  infinite 
periodic  medium  is: 

-  [g(zX) -jp(z.X)]2  E(z.X)  =  0  ( 1 ) 

d 

where  E(z.t,X)=E(z,X)ei',)l  is  an  electric  field  component  (co=  2n|-),  g(z,\)  is  the  gain  (loss) 

constant  and  p(zA)  is  the  propagation  constant.  Both  g  and  (3  are  periodic  in  z  and 
depend  on  \  (the  free  space  wavelength).  (The  dependence  on  X  will  not  always  be 
explicit  in  the  rest  of  this  paper.)  The  usual  assumption  is  that  (g-jJ3)  is  an  even  periodic 
function  of  z.  Most  of  the  review  paper  by  Elachi  [1]  treats  this  case.  Here,  we  will  only 
assume  that  g  and  are  periodic  with  the  same  period  A,  i.eM  that  we  can  write: 

g(z.X)  -  JPfcA) = go  (X) -jfo(X)  +  T  (gq^HPqO-))  exp(j2q-^z)  (2) 

q*0 

(where  g,  and  Pq  are  complex  in  general). 

It  is  well  known  that  the  solutions  of  (1)  are  of  die  Floquet-Bloch  type  [11]: 

E(z)=  A(z)  er* + B(z)  e_r*  (3) 

where  A(z)  and  B(z)  are  periodic  functions  of  z  and  T  is  the  propagation  constant,  which 
is  written  as: 

r=G-jK  (4) 

where  G  is  the  effective  gain  and  K  is  the  effective  propagation  constant.  The  usual 
interpretation  of  (3)  is  to  say  that  A(z)  corresponds  to  a  forward  propagating  mode  and 
B(z)  corresponds  to  a  backward  propagating  mode.  We  shall  see  below  when  this 
interpretation  is  correct.  When  the  Bragg  condition  is  nearly  satisfied  (i.e.,  (3q  close  to 
Kb),  the  waves  reflected  at  each  interface  add  constructively  and  the  reflectivity  becomes 
high.  In  a  passive  medium,  this  gives  a  band  of  frequency  in  which  T  has  a  big  negative 
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real  pan  (as  we  will  see  below),  which  means  that  there  is  no  propagation  in  that  band. 
Hence  the  name  stop-band  is  us  sc.,  which  is  similar  to  the  forbidden  band  in  semiconduc¬ 
tors.  The  Bragg  condition  at  normai  incidence  gives  for  KB: 

Ac  -  *’-3  =  P'r  (5) 

Since  the  band  around  KB  is  die  one  we  are  interested  in,  we  will  look  at  A(z)  and  B(z)  in 
that  band.  Because  A(z)  and  E(z)  correspond  to  independent  solutions  of  (1),  we  can 
treat  them  separately. 

We  start  by  expanding  A(z)  and  B(z)  in  Fourier  series  (which  we  can  do  since  they 
are  periodic): 


A(zA)  =  Ac+  V  Ap(A.)  exp(j2q-^-z) 

n  A 

(6.a) 

q*0 

B(z,A,)asB0+  y  Bq(X)  exp0'2q-7-z) 
n-  -  A 

(6.b) 

9T' 

Then  we  substitute  (2),  (3)  and  (6)  in  (1)  and  by  neglecting  second  order  terms  and  equat¬ 
ing  terms  having  the  same  exponential  factor ,  we  get  after  some  approximations  [6]: 

^(X)  =  G+g0Ij(^)Ao= S+(X)  ^  (7) 

b-W'g^^b°=s-<wb°  (8) 

where: 

6=Kb-Po(X)  5e/f=KB-K 

k.  =  P_p+jg_p  (9) 

(In  equations  (7),  (8),  (9)  and  in  the  following,  p  is  a  positive  integer.)  The  quantities  A^, 
and  B_p  represent  Bragg- scattered  waves  from  A<,  and  B0  respectively.  Their  amplitude 
becomes  significant  only  when  is  near  KB.  We  can  then  neglect  the  other  components, 
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so  that  (3)  becomes: 

E(zX) = Ao  [  1  +  s+(X)  exp(j2KBz)]  er* + B0  [  1  +  s_(X)  exp(-j2KBz)j  e"r*  (10) 


Again,  the  usual  interpretation  is  to  say  that  the  part  in  Aq  is  a  forward  propagating  mode 
and  the  part  in  B0  is  a  backward  propagating  mode. 

Let  us  notice  that  if  (2)  is  an  even  function  of  z,  we  have:  tc=ic_  and  s=s+=s_.  In 
general,  we  have  ls+l  =  ls_l  only  if  f(z)=  g(z)-j(3(z)  is:  (1)  an  even  function  of  z  (k+=  •<_)> 
(2)  an  odd  function  of  z  (<+=-<_).  or  (3)  if  f(z)=  Cgo-j^o)  +(g'-j|3')h(z),  with  h(z)  a  real 
periodic  function  of  z.  This  last  case  gives:  k*=  (g/-jp')cp  and  k_=  (g'-j{3')Cp’  (where  Cp  is 
the  pth  Fourier  coefficient  of  h(z)  and  indicates  the  complex  conjugate  of  Cp),  i.e.: 
I  k*  I  =  l  k_  I .  The  last  case  includes  pure  index  modulation  ^=0)  and  pure  gain  modula¬ 
tion  (|3'=0).  We  will  have  IkJ*Ik_I  (and  thus  i s+ 1  *  I s_ I )  only  if  g(z)  and  (3(z)  have  dif¬ 
ferent  "profiles",  i.e.,  if  f(z)=  (go-j(50)+  [g'h(z)-jP'p(z)]  and  p(z)  is  not  proportional  to  h(z) 
(both  are  real  functions). 

Looking  at  the  scattering  of  and  B_p  into  Ao  and  B0,  respectively,  by  the  same 
method  as  for  (7)  and  (8),  we  get : 


Ao= 


Bn= 


-JK- 


(G-go)+j(5efl-6) 


A+pS  A+p 


-K 


-B__  =  —  B. 


(G-goHj(5efr8) sl“"p 


(11) 

(12) 


Substituting  (7)  into  (11),  or  (8)  into  (12),  we  can  now  obtain  the  characteristic  equation 
(which  also  defines  P): 

P2=(G+-j6eff)2=(g0+j5)2+Kf<_  (13) 


So  that,  by  using  (9)  and  (13)  in  (4): 


T= G+j5efr-jKB = P-jKB 


(14) 
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The  solutions  of  equation  (13)  in  the  case  of  even  functions  are  discussed 

for  example  in  [1],[6]  and  [8].  However,  all  these  analysis  assume  either  a  pure  index 
modulation  or  a  pure  gain  modulation  and  no  material  dispersion  (i.c.,  the  refractive 
index  and  the  gain  are  independent  of  frequency).  This  assumption  is  usually  valid 
because  conventional  edge-emitting  DFB  laser-diodes  have  stop-bands  that  are  2  to  3  nm 
wide,  so  that  index  and  gain  (loss)  variations  can  be  neglected.  But  in  the  case  of  our 
DFB- S ELD  structure  [3],  the  stop-band  is  of  the  order  of  35  nm  and  gain  variations  can¬ 
not  be  neglected  anymore,  especially  when  the  gain  changes  from  negative  to  positive  in 
the  stop-band,  as  we  will  see  below.  Also,  in  the  case  of  our  SELD,  we  have  index  and 
gain  modulation  at  the  same  time  (since  it  is  formed  by  alternating  layers  of  GaAs  and 
GaAlAs  and  gain  exists  only  in  the  GaAs:  see  reference  [2]  for  details). 

Looking  at  (13),  one  question  comes  up:  Which  root  do  we  choose  for  P?  One  of 
these  roots  corresponds  to  the  correct  assignement  of  the  forward  propagating  mode  to 
A(z)  and  the  other  root  makes  A(z)  the  backward  propagating  mode.  This  is  an  impor¬ 
tant  question  since  the  proper  definition  of  the  reflection  and  transmission  coefficients 
depends  on  the  correct  assignement  of  forward  and  backward  propagation  to  the  modes. 
In  the  previous  analysis  [6,8],  P  was  chosen  such  that  G  had  the  same  sign  as  go,  or, 
equivalently,  such  that  5*^  had  the  same  sign  as  5.  But  this  doesn’t  work  in  the  general 
case  and  we  have  to  find  another  way  to  decide  which  root  defines  the  forward  propagat¬ 
ing  mode  as: 

E(z)  =  Ao  (1  +  s„  exp(j2KBz))  e‘ 1  (15) 

Intuitively,  we  see  that  in  the  stop  band  K=KB  (and  5^=0),  so  that  the  direction  of 
the  net  power  flow  can  be  found  by  comparing  the  intensities  of  the  two  components 
(incident  and  backscattered)  in  each  mode.  Looking  at  (15),  it  is  obvious  that  the  direc¬ 
tion  is  reversed  when  I  s+ 1  becomes  greater  than  one. 
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More  rigourously,  if  (15)  is  the  only  non-zero  electric  field  component,  dividing  the 
average  power  flow  per  unit  area  by  the  average  energy  density,  we  obtain  the  energy 
velocity  [6]: 


c  K-lsJ2(2KB-K) 
Ve"  Po  l+ls+l2 


(16) 


where  Po=nciWpo£o<  5  being  the  average  effective  index  and  c  is  the  speed  of  light.  So,  if 
(15)  really  represents  a  forward  wave,  vE  has  to  be  positive.  From  (16),  we  get  then  the 
condition: 


1-t- 1  s.J 2  Kb  _  Kb  _ .  Jirff 
2ls+l2  K  Ka-S^f  Kb 


(17) 


since  in  general  5^  c  KB.  And  this  gives  us: 

I  s+ 1  £  l 


(1 8.a) 


And  similarly,  for  the  term  in  B0  to  be  the  backward  propagating  wave,  we  need  to  have: 


I  s_  I  £  1  (18.b) 

If  we  consider  only  the  cases  when  l$+l  =  ls_l,  there  is  no  contradiction  between  (18.a) 
and  (18.b).  There  is  then  always  a  forward  propagating  mode  and  a  backward  propagat¬ 
ing  mode.  But  if  we  look  at  the  general  case,  it  is  possible  to  have,  for  example  (18.a) 
being  satisfied  while  (18.b)  is  not  satisfied.  To  see  more  clearly  what  happens,  let  us 
define:  P1=G+j5efI  and  p2=-G-j5eff.  Then,  looking  at  (7),(8),(1 1)  and  (12),  we  see  that: 


(Op,* 


1 

(Op, 


(Op,= 


1 

(Op, 


(19) 


where  the  indices  Pi  and  P2  mean  that  the  coefficient  is  computed  using  that  value  for  P. 
Thus,  when  I  s+ 1  =  I  s_  I  =  I  s  I ,  we  have: 
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This  shows  that  in  this  case,  the  choice  of  the  root  is  unambiguous.  Indeed,  if  Pj  gives  a 
Isl  greater  than  one,  P-  will  give  Isl  smaller  than  one  and  vice  versa.  Looking  at  the  physi¬ 
cal  interpretation  of  s+  and  s.  from  (7)  and  (8),  we  see  that  (18)  simply  says  that  the  scat¬ 
tered  wave  amplitude  cannot  be  bigger  than  the  incident  wave  amplitude,  as  we  found 
intuitively  above. 

But  in  the  general  case,  if  1  tc+l  is  different  from  I  tc_l ,  we  will  have  !s+l  different 
from  I  s_  I  and  in  a  certain  frequency  range,  we  could  have  for  example  I  s+ 1  >  1  and 
ls_l<l.  Equations  (19)  then  show  that  this  would  be  true  for  both  roots  of  (13).  If  this  is 
not  an  effect  of  the  approximations  used,  and  there  is  no  reason  to  believe  that  it  is,  what 
does  this  mean?  From  the  discussion  above  about  the  direction  of  the  energy  flow,  the 
interpretation  of  this  result  is  that  both  mode  have  a  net  energy  flow  in  the  same  direc¬ 
tion  (in  this  case,  backward)!  This  does  not  mean  that  there  is  no  energy  going  in  the  for¬ 
ward  direction,  but  only  that  there  is  more  energy  going  backwards,  i.e.  that  the  reflected 
wave  is  greater  than  the  incident  wave  for  one  of  the  modes.  This  is  possible  since  this  is 
an  active  structure. 

The  simplest  structure  that  would  exhibit  this  behavior  is  a  three  layers  periodic 
structure,  for  example:  GaAs/Gai_xAlxAs/  Ga^ALAs,  with  x*y  and  current  injection  to 
have  gain  in  the  GaAs  layer  only.  This  is  because  in  this  case  g(z)  and  (J(z)  would  have  a 
different  "profile".  We  applied  the  theory  to  that  three  layers  structure  and  for  some  gain 
curves,  we  obtained  the  phenomenon  mentioned  above,  i.e.:  in  both  modes  the  net  energy 
flow  is  in  the  same  direction.  We  used  these  results  to  compute  the  reflectivity  and 
transmittivity  of  a  finite  multilayer  structure  [12].  We  then  used  a  transmission  matrix 
method  similar  to  the  one  developed  in  [13]  to  check  the  reflectivity  and  transmittivity 
obtained  with  the  eigenmode  theory.The  agreement  of  the  eigenmode  theory  with  the 
exact  results  was  very  good. 


III.  Two  layers  periodic  structure 

The  case  of  the  two  layers  periodic  structure  is  interesting  for  several  reasons.  First, 
it  is  the  basic  structure  of  our  surface-emitting  laser-diode  (SELD)  [2,3],  which  we  will 
analyze  below.  It  is  also  a  very  common  case  in  applications  such  as  dielectric  mirrors. 
Second,  we  can  get  exact  analytical  solutions  without  too  much  difficulty  and  this  will 
allow  us  to  compare  them  with  those  of  the  approximate  theory  of  section  II  (which  we 
will  call  eigenmode  theory). 

The  structure  used  for  the  calculation  is  shown  in  figure  2.  The  gain  and  propaga¬ 
tion  constants  in  one  period  are  given  by: 

rg1(X)-jp1(X)=r,(X)  if-z,<z<0 

g(zA)-jp(z,X)=|g2(X)_jp2(X)  =  r2(X)  if0<z<Z2  (20) 

and  z!+z2= A,  the  period  of  the  structure.  Let  us  notice  that  if  we  put  the  origin  at  z=— y-, 

the  function  (g(z)-jP(z))  is  even  in  z.  However,  to  simplify  the  calculations,  we  will  use 
the  coordinate  system  of  figure  2  and  correct  later  when  necessary.  We  will  suppose  that 
we  have  a  plane  wave  propagating  in  the  direction  normal  to  the  plane  of  the  layers,  with 
a  linear  polarization  in  the  y  direction,  so  that  equation  (1)  is  the  equation  for  the  y  com¬ 
ponent  of  the  electric  field. 

Several  equivalent  methods  can  be  used  to  solve  equation  (1)  in  this  case.  We  will 
outline  the  one  using  solutions  of  the  Floquet-Bloch  type,  i.e.:  we  start  with  equation  (3). 
(Another  method  is  used  in  [14].)  Since  A(z)  and  B(z)  correspond  to  two  independent 
solutions,  we  can  work  with  them  independently  of  one  another.  We  will  work  out  the 
solutions  for  A(z),  and  B(z)  can  be  obtained  simply  by  replacing  F  by  -I~  everywhere. 
So,  we  substitute  (3)  in  (1)  and  we  get,  for  A(z): 

+  2T-|£  +  (T2  -  (g(z)  -  jp(z))2)  A = 0 

Which  gives  us,  with  (g(z)-j{5(z))  defined  in  (20): 
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j  K\  exp(y+1z)-t- A_i  exp(y_,z) 
A(Z)  =1  A+:  exp(y+2z)  +  A_j  expr^z) 


if  -Zj  <  z  <  0 
if  0<z<Z2 


(21) 


where: 


f  y+i  *-r+r,  7+2=-r+r2 

y_2=-r-r2 

and  A(z)  is  a  periodic  function  of  z  with  period  A=z;+  z2.  Using  the  continuity  of  the 
tangential  components  of  the  electric  and  magnetic  fields  at  z=0  and  between  z=-2l  and 
z=z2  (since  A(z)  is  periodic),  we  get  after  some  algebra: 


>1 

! 

x,1 

j>  >i 

! 

] 

r ,  1 
A+2 

■I'll 

J 

<  1 

i_r2  -nj 

A_2 

h  4 

and: 


SA 


r,e 


.r,i 
-r>er 


e~r'2'  e 
r,*,  _p  „r,zi 


r  1 

A,, 

r,2t  c-r,z« 

r2er*  -r2e~rj! 


'A+2 

A-2 

»  4 

Eliminating  A+2  and  A_2,  we  get: 


(22.a) 


(22.b) 


A  B' 

A+i 

-JH 

C  D 

b  4 

A-i 

W  J 

(23) 


where  p=erA  and: 


_  -riit 


A=e 


B=e 


C=e_r,2< 


Dae"r,z’ 


Il+ 

r2 

f2 

r'J 

r, 

r2l 

r2“ 

ri 

f  r 

1  -L 

r2 

I  r2 

s.  * 

'Tj 

r,  r. 


r,  r. 


sinh(T2Z2)  +  cosh(r2Z2) 

sinh(T2Z2) 

sinh(T2Z2) 

sinh(T 2Z2) + cosh(T 2Z2) 


Which  gives  for  the  eigenvalue  p: 
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Ji  =  erA=y± 


(y)2-l 


(24) 


where  T  is  the  trace  of  the  matrix  and  is  given  by: 


T  =  2cosh(T  i  z  j  )cosh(T  2Z2)  + 


r- 

r2 

— > 

U 

sinh(T  j  z ,  )sinh(T  2Z2> 


In  the  case  of  no  loss  and  no  gain,  this  reduces  to  the  result  obtained  in  [14],  Let  us 
notice  that  the  matrix  in  (23)  is  unimodular  since  AD-BC=1.  This  also  means  that  the 
two  eigenvalues  u.)  and  p2  will  be  inverse  of  one  another 


and 


til  =crA 
ti2=c'rA 


(25) 


We  can  check  that  (24)  satisfies  this  condition.  This  relation  was  expected  for  symmetry 
reasons:  One  eigenvalue  corresponds  to  the  forward  mode  and  the  other  one  to  the  back¬ 
ward  mode.  But  here  again,  we  have  the  same  problem  as  we  had  for  the  approximate 
solution:  Which  one  is  the  forward  mode  and  which  one  is  the  backward  mode?  Since 

we  have  only  erA,  r  is  determined  only  up  to  a  term  of  the  form  (J2m-7-)  (where  m  is  an 

integer).  This  means  that  (T  ±j2m-y)  is  also  a  solution  of  (25),  for  all  m.  The  two  possi¬ 
ble  solutions  for  a  forward  wave  are  obtained  by  choosing  m  such  that  the  wave  vector 
has  a  value  close  to  KB.  We  get  then: 


T =G+j5ei-rjKB  from  Pi  (26.a) 

r=-G-j8efl-jKB  fromp2  (26.b) 

And  we  use  the  same  notation  as  in  section  II.  We  see  that  this  is  totally  equivalent  to  the 
choice  of  the  root  of  equation  (13). 

Hence,  we  will  use  the  same  criterion  as  in  section  D:  the  direction  of  the  net  power 
flow.  To  do  that,  we  first  have  to  compute  the  eigenvectors  corresponding  to  the  eigen¬ 
values.  We  start  by  rewriting  (21)  as: 
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J  A+,(exp(Y+,z>+s,exp(y_,z))  if-z^zcO 
A(z j  1  (S3  exp(Y+Jz)  +  s4  exp(Y_2z))  if  0  <  z  <  z; 


I 


(27) 


Then,  from  (23),  we  get: 


,  A-i  n-A  C 
Sl"  A+i  ”  B  %-D 


(23) 


And,  with  (22.a): 


S3  =  Ti  =  ^(ri+r^'(r2"r>)) 

A+j  2r2 

S4 =  ~r~  =  T7r'CTi-Tr,^i'T2+r  1)) 

a+i  ir2 


A_2 

Notice  that  if  we  define:  s2-- — = 

A+2 


S4 

— ,  we  get: 

s3 


st+a 

S2~  1+asj 


r*-r» 

w“h: 


(29) 


This  is  a  biiinear  transformation  in  the  complex  plane,  and  if  a*0  (i.e.:  ri*r2),  it  maps  the 
disk  I  si  I  <1  onto  another  disk,  but  with  a  part  that  has  I S2 1  >1.  We  can  thus  not  rely  on  si 
or  s2  alone  to  get  the  direction  of  the  power  flow,  in  the  way  we  looked  at  s  in  section  II. 

To  find  the  direction  of  the  power  flow,  we  will  look  at  the  Poynting  vector’s  direc¬ 
tion.  But,  because  we  have  gain  and/or  loss  in  the  structure,  the  Poynting  vector  is  not 
constant  along  z.  Therefore,  we  will  take  an  average  over  one  spatial  period  (of  the  time 
averaged  vector): 

<S>  =  —  f  i-Re(£xH*)dz  (30) 

A  '  2 


where  we  take:  Ey=A(z)exp(Tz)  and  A(z)  is  defined  in  (27).  An  explicit  formula  is 
derived  in  Appendix  A  and  it  shows  that  the  only  non-zero  component  of  the  Poynting 
vector  is  <S^>.  The  of  <S*>  will  then  give  us  the  direction  of  the  power  flow:  posi¬ 
tive  means  forward  propagation  and  negative  means  backward  propagation. 
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Another  way  to  find  the  direction  of  the  power  flow  is  to  find  the  s  coefficient 
defined  in  section  H  To  find  s,  we  compute  the  first  two  coefficients  of  the  Fourier  series 
expansion  of  A(z)  and  get  s  as  their  ratio.  We  can  then  look  at  Isl  as  in  section  II.  How- 

ever,  we  have  to  be  careful  about  one  thing:  the  point  of  even  symmetry  is  at  z=-—  and 

that  is  where  we  have  to  put  the  origin  when  we  compute  s.  If  we  do  not  do  that,  we  will 
have  a  phase  difference  between  the  s  coefficient  computed  here  and  the  one  defined  in 
section  n.  Also,  we  would  not  have  s+=s_=s.  Taking  that  into  account,  we  get : 


s=e“jK*2' 


j  A'(z)  dz  /  j  A(z)  dz 


“I I 


(31) 


where  the  exponential  factor  is  due  to  the  change  of  origin,  A(z)  is  given  by  (27)  and 
A'(z)  also,  but  with  the  y^  replaced  by  y^'  (i=l,2): 


Y±,'=Y±i-j2KB 


Explicitely: 


jA(z)dz= 

-t\ 


Y+i 


+  s, 


Y-i 


+  S3 


Y+2 


+  S4 


eT~aZ“-l 

Y-2 


(32) 


And  the  same  for  A'(z),  but  with  the  y^'. 

A  comparison  of  the  results  obtained  with  lskl  and  with  the  average  of  the  Poynting 
vector  on  the  examples  of  the  next  section  shows  only  negligible  differences. 


IV.  Examples  and  discussion 

To  illustrate  the  theory  outlined  above,  we  will  use  the  two- layer  structure  of  our 
surface-emitting  laser-diode  [2].  The  layer  structure  is  shown  in  figure  2  and  the  parame¬ 
ters  are  given  in  table  1.  Basically,  we  have  alternating  layers  of  GaAs  and  GaojAlo  jAs. 
The  thickness  of  each  layer  is  chosen  to  be  where  is  the  wavelength  at  which 
we  want  maximum  reflectivity  and  n;  is  the  index  of  refraction  of  the  layer  (at  Xo).  The 


GaAs  layer  will  be  layer  nur.  r  1  and  the  GaAlAs  layer  will  be  layer  number  2.  The 
loss  in  the  GaAlAs  layer  is  taken  as  independent  of  frequency  and  equal  to  10  cm-1 
(intensity  loss).  The  gain  in  the  GaAs  layer  depends  on  the  carrier  density.  A  semi- 
empirical  model  (based  on  absorption  curves  in  [15])  was  developed  to  get  analytical 
expressions  for  the  gain  as  a  function  of  carrier  density,  temperature  and  wavelength. 
Figure  3  shows  typical  gain  curves  for  three  different  carrier  densities  at  300  K.  These 
will  be  used  later  in  the  examples. 

Notice  that  the  coupling  constant  for  our  structure  (k=  4500  cm-1)  is  two  orders  of 
magnitude  larger  than  for  edge-emitting  semiconductor  DF3  lasers.  Also,  since  only  the 
GaAs  layer  will  be  pumped,  we  have  not  only  a  periodic  variation  of  the  index  of  refrac¬ 
tion,  but  also  of  the  gain  constant.  This  structure  has  even  symmetry  and  thus  s=s+=s_. 

Figure  4  shows  the  effect  on  G  and  Seff  of  the  introduction  of  gain  and  then  gain 
modulation  in  this  structure.  In  fig.  4(a),  there  is  no  loss  or  gain  and  we  see  that  G  is  zero 
everywhere,  except  in  the  stop-band,  where  it  is  negative  (indicating  total  reflection,  not 
absorption).  In  fig.  4(b),  a  constant  (independent  of  frequency  and  of  z)  gain  of  2000 
cm-1  is  introduced  in  the  structure.  And  in  fig.  4(c),  the  periodicity  of  the  gain  is  taken 
into  account  (i.e.:  we  have  gain  only  in  the  GaAs  layer,  but  with  the  same  average  gain). 
We  see  that  the  effect  of  gain  periodicity  cannot  be  ignored. 

'  Figure  5  shows  G  and  5^  for  different  carrier  densities.  In  this  case,  the  gain  is  not 
constant  but  depends  on  the  wavelength.  The  intensity  gain  in  the  GaAs  layers  used  for 
these  examples  is  shown  in  figure  3.  These  calculations  were  done  using  the  approximate 
(eigenmode)  method  of  section  II  and  the  exact  method  of  section  HL  The  results  of  the 
two  methods  for  G  and  5^  are  undistinguishabie  at  the  scale  of  this  figure.  The  most 
surprising  feature  of  these  results  is  the  sudden  jump  of  G  from  a  negative  to  a  positive 
value  when  N=1.0xl0I8cnT3  (fig.  5(b)).  Notice  that  this  jump  occurs  at  the  wavelength 
where  the  average  gain  (go)  changes  from  negative  to  positive. 
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But  this  is  not  as  surprising  as  it  seems.  If  we  look  again  at  fig.  4(a)  and  4(b),  we 
see  that  when  the  gain  is  changed  from  negative  to  positive,  G  (in  the  stop-band)  changes 
from  a  large  negative  value  to  a  large  positive  value  without  going  to  zero.  What  hap¬ 
pens  here,  in  fig.  5(b),  is  that  the  gain  changes  sign  in  the  stop-band  and  thus  G  changes 
sign  also.  Mathematically,  Isl  becomes  greater  than  one  and  thus  we  have  to  switch  to  the 
other  solution  of  (24)  (or  equivalently  of  (13)),  which  implies  a  sign  change  of  G.  Figure 
6  shows  a  plot  of  Isi  as  a  function  of  wavelength  for  the  two  possible  solutions  ((26.a)  and 
(26.b)).  Since  lsl>l  means  that  the  direction  of  the  power  flow  is  reversed,  we  see  that 
the  eigenmodes  exchange  their  roles  at  the  point  where  lsl=l:  The  mode  that  was  pro¬ 
pagating  in  the  forward  direction  now  becomes  the  backward  propagating  mode  and 
vice-versa. 

Physically,  as  we  can  see  from  (16),  the  velocity  of  energy  propagation  is  very  small 
in  the  stop-band  (Isl  is  close  to  one),  because  of  the  multiple  reflections.  Thus  a  small 
loss  or  a  small  gain  will  be  amplified  a  lot  by  the  multiple  reflections.  In  fact,  we  have 
[6]: 


where  v  is  the  velocity  of  light  in  the  medium  (without  the  periodicity)  and  vE  is  the 
energy  velocity  defined  in  (16).  But  this  jump  is  only  a  mathematical  artifact  due  to  the 
method  used  here.  The  rate  of  change  of  the  energy  density  (which  is  equal  to  gov=  GvE, 
see  [6])  is  a  continuous  function  of  wavelength.  If  we  compute  the  reflectivity  or 
transmitivity  of  a  finite  multilayer  with  these  formulas,  we  get  a  perfectly  continuous 
curve  (as  we  will  see  below). 

Figures  5  and  6  also  show  that  for  N=  0.2xi018cm“3  and  N=  2.0xl018cm-3,  this 
phenomenon  does  not  occur.  For  N=  0.2xl0I8cm"3,  the  gain  is  always  negative  and  we 
don’t  have  any  problem.  For  N=  2.0x1 018cm"3,  the  sign  change  of  the  gain  (go)  occurs 
outside  of  the  stop-band  and  we  see  that  G  changes  sign  also,  but  without  any 


discontinuity.  Figures  5(a)  and  5(c)  show  that,  for  these  cases,  there  is  no  exchange  of 
roles  between  Lie  two  modes. 


In  figure  7,  we  compare  the  ted  and  imaginary  part  of  s,  as  computed  by  the 
approximate  method  (eq.«7i)  -ind  by  the  exact  method  (eq.(31)).  We  see  that  the  differ¬ 
ence  is  small  and  that  for  N=  1.0  *lCI8cnT',  the  discontinuity  occurs  at  exactly  the  same 
wave'ength.  From  this  and  the  fact  that  that  G  ;tnd  8^  are  correct  within  a  few  percents 
with  ihe  approximate  method,  we  can  conclude  that  the  approximate  method  is  a  good 
approximation.  We  are  interested  in  the  s  coefficient  because  it  will  be  used  in  the  next 
section  tc  compute  reflection  coefficients  and  the  emission  spectrum  of  our  DFB-SELD. 


V.  Reflection  and  transmission  coefficients 

We  saw  above  that  for  a  periodic  one-dimensional  medium  with  even  symmetry,  a 
gocd  approximation  of  the  electric  field  is  given  by: 

•  E(z) = Ao(l-!-€f)ert+  Bo(I+sb)e“r*  (34.a) 

where  A<,  and  30  are  the  amplitude  coefficients  for  the  forward  and  backward  propagating 
modes  respectively  and: 


Sf=sexp(j2KBz) 

sb=sexp(-j2KBz) 


(34.b) 


T,  Kb  and  s  have  been  defined  in  equations  (4),  (5)  and  (7),  respectively. 

From  these  definitions,  we  can  immediately  find  the  forward  propagation  factor 
from  z,  to  z2  (with  Zjcz^ffi],  that  is  the  factor  Df  such  that  E(Z2>=Df  E(z,) : 


1  +t  exp(j2KBz2)  a_2+sn  a 
l+sexp(j2KBz,)  ”  1+Sf| 


(35.a) 


where  L^z2-z].  Similarly,  for  the  backward  propagating  factor  Dj,  (from  z2  to  z,),  such 


that  E(z,)=DbE(z2),  we  get: 


l-j-scxpf-^KflZ!)  a_i±f*L  rx 
1  +  s  exp(-j2KBzj)  “  1  +  sb2 


I 


(35.b) 
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One  important  point  when  using  these  formulas  is  that  the  origin  of  the  z  axis  has  to  be  at 
a  point  of  even  symmetry.  Otherwise,  we  would  have  made  a  mistake  because  the  s 
coefficient  for  the  forward  and  backward  modes  would  have  a  different  phase  (see  the 
definitions  (7)  and  (8)). 

A.  One  interface 

Let  us  now  look  at  what  happens  at  the  interface  between  a  periodic  region  and  a 
uniform  region  (Fig. 8).  Expression  (34)  for  the  electric  field  in  a  periodic  region  is  only 
an  approximation,  but,  as  we  saw  above,  it  is  a  good  approximation.  We  will  thus  sup¬ 
pose  in  the  rest  of  this  paper  that  the  electric  field  in  a  periodic  region  is  of  the  form 
given  by  (34).  In  the  uniform  region,  the  electric  field  will  be  given  by: 

Eu(z)  =  A'0e~jr,(*'I,) + B,0e)T,(*~z')  (36) 

where  ru  is  the  complex  propagation  constant  in  the  uniform  region  and  A'0  and  B'0  arc 
the  amplitudes  of  the  forward  and  backward  propagating  modes  respectively,  zj  is  the 
position  of  the  interface,  which  we  can  subtract  from  z  in  (36)  to  simplify  the  calcula¬ 
tions,  since,  for  a  uniform  region,  the  choice  of  the  origin  is  arbitrary. 

Figure  8  shows  the  meaning  of  the  different  reflection  and  transmission  coefficients 
we  are  going  to  compute.  To  find  these  coefficients,  we  express  the  continuity  of  the 
tangential  electric  and  magnetic  fields  at  z=zj  and  we  get: 

(1  +  sfl)A0eri‘  +  (1  +  sbl)Boe-rz'  =  A'0 + B'0 

(Sf!  -  1)Kb  Aoerit +( 1  -  s^KsBoe-^ = -&A'0 + puB'0  (37) 

where  the  index  1  for  Sf  and  sB  indicates  that  they  are  computed  at  z=z,  and  where  we 
used  T  =  Kb  (we  are  working  near  the  Bragg  wavelength)  and  Tu  =  pu.  Now,  to  find,  for 
example,  R!(  we  have  to  take  the  ratio  of  the  backward  and  forward  propagating  electric 
fields  at  the  interface  in  the  periodic  region: 

R  Ep,b«*(zi)  B0(l+sbl)e~r*' 

Ep/orw(z  t )  Aq(  1  +  Sfj  W1' 


(38) 
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when  there  is  no  node  incident  from  the  right  (B'0=0).  We  solve  then  (37)  to  find  Bq/Aq 
and  plug  the  result  in  (38).  The  other  coefficients  are  found  in  a  similar  way.  We  get: 

„  r-sfl  l+sfc]  ^  sbl-r 

K>  = - -  K?  s  ■" 

1-rSbi  1+Sfi  1  —  rsbi 


T, 


=  (l+r) 


1  ~  sri^bi 

(i  —  rSbj)(l  +Sfi) 


T2  =  (l-r) 


1-rSb, 


(39) 


where  r  is  the  reflection  coefficient  due  to  the  change  in  propagation  constant: 

Kb-Pu  n-nu 
KB  +  PU  n  +  rtu 


(40) 


where  n  is  the  effective  index  of  refraction  in  the  periodic  region  at  the  Bragg  wavelength 
(n  =  Kb/Icq  and  ko=2jt/\).  and  n^  is  the  index  of  refraction  of  the  uniform  region  at  the 
same  wavelength. 

We  see  that  in  the  limit  of  no  periodicity  (s  going  to  zero),  we  recover  the  usual 
expressions  for  the  reflection  and  transmission  coefficients.  And  in  the  limit  r  going  to 
zero,  we  get  the  same  results  as  in  [6].  Now,  if  the  periodic  region  is  on  the  right  instead 
of  on  the  left  as  in  figure  8,  we  just  have  to  exchange  Sf  with  Sb,  Rj  with  R2  and  Tj  with  T2 
in  (39)  to  get  the  formulas  for  this  case.  If  we  call  R3  the  new  R,  we  obtain  in  this  way, 
R4  the  new  R2  and  so  on  and  if  we  replace  Zj  by  z'2,  we  find: 


sW 

I-i's'h 


I'-S'g  1+s'g 
l“f  s' &  1  +  S^ 


T3 


=d-0 


l+gc 

l-rYn 


T4=(l+0 


1  —  s'jgS'p 

(i-r's'{2)(i+s'b2) 


(41) 


I 


B.  Phase  shifter 


In  the  following,  we  will  also  need  expressions  for  global  reflection  (R\  R")  and 
transmission  (T,  T")  coefficients  of  a  phase  shifter.  By  phase  shifter,  we  mean,  in  this 
paper,  a  uniform  region  between  two  periodic  regions  (see  figure  9).  The  Rj’s  and  T,’s 
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(i=l,2,3,4)  are  the  coefficients  given  by  (39)  and  (41).  Notice  that  we  are  using  two  dif¬ 
ferent  coordinate  systems.  This  is  necessary  because,  as  mentioned  earlier,  we  have  to 
put  the  origin  at  a  point  of  even  symmetry  in  a  periodic  region.  But,  in  general,  it  is 
impossible  to  have  the  origin  (cr  its  translation  by  a  multiple  of  the  period  A)  at  a  point  of 
even  symmetry  in  the  two  periodic  regions  at  the  same  time.  Hence  the  use  of  z  for  the 
left  region  and  z  for  the  right  region. 

The  propagation  factor  for  a  uniform  section  of  lenght  and  with  propagation  con¬ 
stant  Tu=pu+j  go  is: 

Pu=exp(-jruLu)  (42) 


We  will  suppose  here  that  the  reflection  coefficient  r  is  equal  to  zero  (i.e.,  n=nu,  which  is 
approximately  the  case  in  our  laser  structure).  We  get  then,  using  the  method  of  multiple 
reflections: 


Pu2s'a-sn 

1  +  sbl 

l-P^s'n 

1+Sfi 

Pu(l-SflSt,l) 

l+s'G 

1  _Pusbls,f2 

1+Sfl 

Pu  sbl  -  S'b2 

l+s'a 

1-PuSbiS'n 

1  +S'b2 

P„(  1  -  S' (2^ b2)  1  +Sbi 

1 -Pu2sbIs'f2  1+s'b2 


(43. a) 


(43.b) 


(43  .c) 


(43  .d) 


The  prime  on  the  s  indicates  that  the  s  coefficient  for  the  right  region  can  be  different 
from  the  one  for  the  left  region. 


VI.  Reflectivity 

Figure  10  and  table  1  show  the  layer  structure  used  for  our  surface  emitting  laser 
diode  (SELD).  From  the  top,  we  have  20  pairs  of  layers,  consisting  of  one  GaAs  layer 
and  one  GaAlAs  layer  (each  layer  is  one  quarter  of  a  wavelength  thick),  then  we  have  a 
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quarter-wavelength  GaAs  layer  acting  as  phase-shifter  and  finally  60  more  pairs  of 
GaAs/GaAlAs  above  the  GaAs  substrate  (see  [2,3]).  This  structure  was  chosen  because 
it  allows  lasing  in  the  middle  of  the  stop-band  [16],  and  thereby  avoids  the  degeneracy 
problem  of  usual  DFB  lasers.  The  number  of  pairs  at  the  top  is  smaller  so  that  most  of 
the  power  output  is  through  the  surface  and  not  in  the  substrate.  The  structure  was  grown 
by  metallo-organic  chemical  vapor  deposition  (MOCVD)  for  the  devices  we  will  present 
here,  but  we  have  also  several  wafers  grown  by  molecular  beam  epitaxy  (MBE). 

In  this  calculation,  we  will  consider  only  plane  waves  at  normal  incidence,  coming 
from  the  air.  Figure  11  shows  the  model  used  to  compute  the  reflectivity.  L(  is  the 
length  of  the  20  pairs  section,  U,  is  the  length  cf  lie  phase  shifter  and  L3  is  the  length  of 
the  60  pairs  section.  The  Rj’s  and  the  T,’s  we  need  are  given  by  the  formulas  derived  in 
the  previous  section  and  are  summarized  in  Appendix  B. 

Using  the  method  of  multiple  reflections,  we  find  the  total  power  reflectivity  to  be 
given  by: 

RP=  I  Rt  1 2  (44.a) 

and 


Rt = Rj  + 


Reff=R3  + 


Tl^RrffPiP; 
1  -R2Re£fPlP2 
T3T4R5P3P4 
1  -R4R5P3P4 


(44.b) 


where  the  P^s  are  propagation  factors  of  the  type  given  by  (35)  (Df  and  Db)  and  R,^  is  an 
effective  reflection  coefficient  for  the  phase-shifter  and  everything  that  is  to  its  right  (in 
fig.  II).  If  we  look  at  the  formulas  we  obtained,  we  see  that  a  lot  of  factors  of  the  form 
(1  +Sf)  or  (1  +Sb)  can  be  cancelled  out.  This  can  be  done  by  omitting  the  factors  between 
square  brackets  in  Appendix  B. 

Before  we  can  compute  the  power  reflectivity  of  our  structure,  we  need  to  define 
two  more  characteristics:  the  index  of  refraction  and  the  absorption  as  functions  of 
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wavelength.  We  will  assume  that  the  index  of  refraction  is  constant  and  that  the  absorp¬ 
tion  coefficient  in  the  GaAlAs  layer  is  constant  and  equal  to  10  cm-1  (intensity 
coefficient).  For  the  GaAs  layer,  we  will  use: 

r  '\~l 

o(E)=  ±  +  ±  [cm'1] 

P  Y 

v.  ^ 

|5=exp(4.85E+ 1.9552) 

y=exp(118.84E- 159  9106)  (45) 

where  E  is  the  photon  energy  in  eV.  This  formula  was  obtained  by  fitting  to  a  measured 
absorption  curve  for  small  p-tvpe  doping  given  in  [15]. 

Now  we  can  compute  the  power  reflectivity  of  our  structure.  Figure  12(a)  gives  the 
result  of  the  calculation  for  the  structure  as  defined  by  figure  10  and  Table  1.  In  Figure 
12(b),  we  see  a  typical  experimental  result  The  measurements  were  made  using  a 
Perkin-Elmer  Lambda  9  spectrometer  with  the  Absolute  Reflectance  Accessory.  In  this 
apparatus,  the  beam  normally  reflects  twice  from  the  sample  at  two  different  positions, 
with  an  angle  of  8  degrees.  Because  of  the  small  size  of  the  sample  and  of  its  non- 
uniformity,  an  aluminium  miiTor  was  substituted  at  the  second  reflection  point.  We  did  a 
first  measurement  with  the  sample,  then  replaced  it  with  an  aluminium  reference  mirror. 
This  allowed  us  to  compute  the  absolute  reflectivity,  knowing  the  reflectivity  of  the 
aluminium.  The  fact  that  we  were  not  at  normal  incidence  (hut  8  degrees  off)  slightly 
increases  the  apparent  thickness  of  the  layers  Oess  than  1%),  but  can  otherwise  be 
neglected. 

We  can  see  that  the  experimental  result  is  similar  to  the  theoretical  result  of  fig. 
12(a),  but  that  there  are  some  differences  also.  The  main  difference  is  the  double  dip  in 
the  middle  of  the  stop-band  for  the  experimental  result  This  can  be  explained  by  looking 
at  what  happens  when  the  thickness  of  each  layer  changes  by  one  lattice  constant  (a=5.65 
Angstroms).  This  gives  us  a  change  of  tickness  of  about  0.9  %,  which,  for  an  original 
Bragg  wavelength  of  890  nm  (in  vacuum),  corresponds  to  a  change  of  about  8  nm.  This 
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corrcsponds  approximative^  to  the  spacing  of  the  two  dips.  What  happens  is  that  the 
grown  structure  is  not  uniform  across  the  wafer  and  the  spot  size  tor  the  reflectivity 
measurement  is  several  miliineters.  Tne  measured  reflectivity  is  thus  an  average  of  two 
different  thic<tnesses. 

We  varied  the  parameters  slightly  to  try  to  reproduce  the  experimental  curve.  Our 
best  result  is  the  theoretical  curve  in  Fig.  12(b)  (but  it  is  not  an  optimal  fit).  It  was 
obtained  by  a  weighted  average  between  a  Bragg  wavelength  (Ao)  of  889  nm  (70  %)  and 
a  Xq=  897  nm  (30  <£),  with  the  thickness  of  the  top  GaAs  layer  (the  one  in  contact  with 
air)  reduced  to  37  nm  (instead  of  o7  nm).  (This  reduction  in  thickness  could  happen  dur¬ 
ing  wafer  processing  and  handlmg.) 

We  did  not  try  to  improve  the  agreement  between  theory  and  experiment  because  it 

depends  on  too  many  unknown  parameters,  such  as  the  exact  dependence  of  the  loss  on 

wavelength,  and  because  the  theory  used  here  is  only  an  approximation.  We  also 

• 

neglected  the  variation  of  the  index  of  refraction  with  wavelength.  However,  we  get 
correctly  the  main  characteristics  of  the  curves  and  this  gives  us  confidence  for  applying 
the  theory  to  the  calculation  of  the  emission  spectrum  of  the  SELD  (in  section  Vm). 

VII.  The  surface-emitting  laser  diode 

As  described  previously  in  [2]  and  [3],  the  surface  emitting  laser  diodes  (SELD)  are 
fabricated  by  starting  from  wafers  having  the  layer  structure  described  above  (Fig.  10  and 
Table  1).  Mesas  with  dimensions  ranging  from  2  x  8  to  3  x  17  pm2  at  the  top  surface  are 
formed  by  wet  etching.  Then  a  n-type  Gao.6Alo.4As  cladding  layer  is  grown  around  the 
mesas  by  selective  liquid  phase  epitaxy  (LPE).  Finally,  a  lateral  pn  junction  is  formed  by 
selective  zinc  diffusion  through  an  opening  in  a  silicon  nitride  film,  at  a  distance  of  3  to  4 
pm  away  from  the  edge  of  the  mesa.  Figure  13(a)  shows  the  structure  of  the  laser  and 
Fig.  13(b)  is  a  scanning  electron  microscope  (SEM)  picture  of  a  laser.  The  zinc  diffusion 
region  is  brigther  in  this  picture.  The  total  thickness  of  the  structure  is  about  10  pm. 
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Since  the  GaAs/GaAlAs  multilayer  is  completely  surrounded  by  GaAlAs,  carrier 
confinement  in  this  structure  should  be  comparable  with  that  of  buried  heterostructures. 
Carriers  are  injected  predominantly  in  the  GaAs  layers  and  therefore  we  will  assume  in 
the  model  that  only  the  GaAs  layers  are  pumped.  There  is  however  the  possibility  of 
current  leakage  around  the  mesa.  We  believe  this  happens  at  high  currents. 

Figure  14  shows  the  light  output  versus  dc  current  characteristics  (L-I  curve)  of  a 
SELD  at  room  temperature.  We  get  what  looks  like  typical  threshold  currents  in  the 
range  from  2  to  10  mA  and  differential  quantum  efficiencies  from  0.3  %  to  1.36  %. 
Powers  of  0.5  mW  for  a  current  of  50  mA  have  been  observed  for  CW  operation  at  room 
temperature.  More  recently,  we  obtained  quantum  efficiencies  up  to  5  %  for  some  dev¬ 
ices.  Unfortunately,  they  were  accidentely  destroyed  before  their  spectrum  could  be 
measured. 

Observation  of  the  near-field  pattern  shows  [2,3]  that  the  light-emitting  region  is 
confined  within  the  rectangular  top  surface  of  the  multilayer  (at  least  when  the  current  is 
not  too  high).  The  far-field  pattern  is  circular  [2,3]  and  its  beam  divergence  angle,  for  a 
current  of  20  mA,  is  7  to  8  degrees.  However  the  circular  emission  pattern  is  surprising 
for  a  device  with  a  rectangular  emission  surface.  The  7  to  8  degrees  angle  corresponds 
roughly  to  a  width  of  about  3  pm  for  the  rectangular  opening,  but  the  length  of  about  15 
pm  should  give  a  much  narrower  beam  in  that  direction  (about  3  degrees),  if  it  were  a 
diffraction  limited  beam.  Obviously,  the  beam  is  not  diffraction  limited  in  that  direction. 
A  reason  for  that  could  be  that  there  is  no  coherence  between  points  more  than  3  pm 
apart. 

Figure  15  shows  the  spectrum  of  the  SELD  at  different  currents,  for  CW  operation 
at  room  temperature.  We  see  that  there  is  a  sharp  emission  line  at  884  nm  which 
becomes  obvious  for  a  drive  current  of  1.6  mA.  This  is  consistent  with  the  L-I  curve, 
which  shows  what  looks  like  a  2  mA  threshold  current.  The  half-width  is  rather  large 
(1.5  to  3  nm)  compared  to  conventional  edge-emitting  laser  diodes.  The  total  width  of 
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the  emission  spectrum  is  also  large  (about  50  nm),  compared  to  edge-emitting  lasers. 

In  view  of  these  results,  one  might  ask  whether  these  devices  are  really  lasing?  We 
think  they  are.  There  may  be  several  reasons  for  the  low  quantum  efficiency:  leakage 
current  around  the  mesa,  non-radiative  recombination  at  interfacial  defects  if  the  L?E 
was  not  very  good  and  a  big  spontaneous  emission  component.  We  explain  below  how 
the  shape  of  the  mesa  can  lead  to  a  lot  of  parasitic  spontaneous  emission  which  lower  the 
quantum  efficiency  and  can  also  affect  the  spectrum. 

VIII.  Emission  spectrum  of  the  SELD 

In  this  section,  we  develop  a  modei  (figure  16)  to  reproduce  and  explain  the  emis¬ 
sion  spectrum,  at  least  below  threshold  It  is  very  similar  to  the  one  used  to  compute  the 
reflectivity  and  we  make  the  same  assumptions,  except  that  we  will  now  have  gain  in  the 
pumped  section  (which  is  taken  as  being  the  zinc  diffused  region).  In  addition,  we 
assume  that  this  pumped  section  has  a  uniform  carrier  density  (in  the  GaAs  pan  of  the 
layers)  and  that  the  unpumped  region  is  uniform  and  lossy.  These  assumptions  are  rough 
approximations,  but  they  allow  us  a  relatively  easy  calculation  using  what  we  did  in  the 
previous  sections. 

The  reflection  and  transmission  coefficients  are  given  in  Appendix  B.  Eb  E4  and  Ej, 
Ej  are  the  electric  fields  amplitudes  of  the  forward  and  backward  propagating  modes 
respectively,  at  the  two  interfaces  of  the  phase  shifter  region  (Lj).  E^  is  the  amplitude  of 
the  electric  field  of  the  light  output  The  Sj’s  are  the  equivalent  sources  for  the  electric 
fields  at  the  interfaces  due  to  the  spontaneous  emission  in  sections  L]  and  L3  (we  neglect 
the  spontaneous  emission  in  the  phase  shifter). 

If  we  express  the  electric  fields  as  functions  of  one  another  and  of  the  S/s,  we  get 
the  following  system  of  equations: 

Ei  =R2PiP2E2+R2PiS1  +  S2 
E2  =  R3E|  +T4E5 

—  RSP3P4E4  +  R  jP4S4  +  S3 
E4  =  R4E3  +T3Ej 
E.^=T1PiE1+TJS< 


(46) 
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If  we  eliminate  the  Ej’s,  we  find: 

EoUt=^;T2P,P2^  +  T^S ,  +  (T2P2  £  )S2 + (T2P2  ^)S3  +  (R5T2P2P4  ^-)S4  (47) 


where: 


N  =  R3 +RjP3p4(T3T4  -R3R4) 

D=(l-R2R3PiP2)(l-R4R5P3P4)-R2R5T3T4P1P2P3P4 

But  what  we  are  interested  in  is  the  output  power.  To  convert  the  relation  (47)  between 
the  electric  fields  to  a  relation  between  the  power  in  the  modes,  we  have  to  remember 
that  the  S^s  are  really  of  the  form  given  by  (34),  i.e.: 

Si = S  i0(  l+s  ex*±  2jKBZi))  e^ 

while  we  have: 


E^z)=Eoe  and  E^g — E^^Z]) = E^ 

Then,  with  the  same  cancellations  of  factors  as  for  the  reflectivities,  we  get: 

I Eout  1 2  =  uw^p^+t^2  is10i2+  it2p2£i2  is20i2+  IT2P2^-I2  I S30 1 2 

+  IR5T2P2P4^-12IS40I2  =  £cj  isj0l2  (48) 

u  j-l 

where  we  have  taken  into  account  the  fact  that  the  S^s  are  not  coherent  with  one  another. 
This  formula  is  valid  with  the  reflection  and  transmission  coefficients  as  given  in  Appen¬ 
dix  B  when  we  omit  the  factors  between  square  brackets.  Now,  by  computing  the  Poynt- 
ing  vector  for  a  field  of  the  type  of  (34),  we  find  a  relation  between  the  power  density  p; 
and  !  Sj0 1 2: 

Pj  =  2^^  CK“  ‘ S  ‘  2(2Kb“K)1  1  Sj° 1 2 


And  for  the  output  power: 


ko 

P”“"  20)40 


lEol2 
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This  gives  us  then  the  ouput  power  spectral  density: 


PoUl =  . 


ko 


K-IsI^Kb-K)  \~[  J‘Jj 


f  4  1 

i  sciPjr 


(49) 


v/here  the  Cj’s  were  defined  in  (48)  and  the  p^s  (defined  below)  are  the  spectral  power 
densities  due  to  the  spontaneous  emission. 

We  see  immediatly  from  (48)  and  (49)  that  we  will  have  lasing  when  D=0.  This 
condition  will  give  us  the  wavelengths  of  the  modes  and  the  corresponding  gain  thres¬ 
holds.  (Do  not  forget  that  the  reflection  and  transmission  coefficients  are  wavelength 
dependent  and  that  the  propagation  factors  depend  on  the  wavelength  and  on  the  gain.) 

If  we  solve  the  equation  D=0,  we  get  several  solutions,  each  one  corresponding  to  a 
mode.  If  we  use  Xo=  890  nm  and  a  phase  shifter  thickness  reduced  to  0.4  A  (instead  of 
0.5  A),  we  will  see  below  that  we  get  the  best  spectrum  adjustment.  This  gives  us  a  mode 
at  884  nm  with  a  threshold  intensity  gain  of  about  44  cm-1  (for  the  material  gain  of 
GaAs).  With  the  model  for  the  gain  developed  below,  this  corresponds,  at  that 
wavelength  to  a  threshold  carrier  density  of  about  1.6xl018  cm"\  We  can  now  estimate 
the  threshold  current  with: 


tldN* 

where  t=  2  pm  and  1=  15  pm  are  the  mesa  dimensions  and  d=  3  pm  is  the  active  layers’ 
thickness  (without  the  GaAlAs  layers)  and  t=  4  ns  is  the  spontaneous  lifetime  (q  is  the 
electron  charge).  This  gives  us  =  7.76  mA,  which  is  in  the  range  of  observed  threshold 
currents.  The  other  modes,  which  will  give  humps  in  the  spectrum,  have  much  higher 
threshold  gains  (several  thousands  cm-1  at  least). 

Let  us  now  find  expressions  for  the  pfs  in  terms  of  the  spontaneous  emission  spec¬ 
tral  density  rate  and  the  device  parameters.  To  simplify  the  calculations,  we  will  suppose 
a  uniform  spontaneous  emission  in  the  whole  section  (which  is  consistent  with  the  uni¬ 
form  carrier  density)  and  we  correct  then  for  the  fact  that  we  have  emission  only  in  the 
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GaAs  layers  with  a  factor  f.  This  factor  f  is  equal  to  the  ratio  of  the  thickness  of  one 
GaAs  layer  to  one  period.  Now,  to  get  the  p^s,  we  sum  incoherently  the  amplified  spon¬ 
taneous  emission  in  each  section.  For  example,  we  get  for  p^ 

L“ 

Pl  =  fT-Er,p(E)Y  Jlerzl2dz 


where  y  is  the  fraction  of  the  spontaneous  emission  coupling  to  the  cavity  mode,  E  is  the 
photon  energy  (in  Joules)  and  r^,(E)  is  the  spontaneous  emission  spectral  density  rate  per 
unit  volume.  The  cavity  cross-section,  txl,  is  taken  as  constant  The  factor  one  half 
comes  from  the  fact  that  only  one  half  of  the  spontaneous  emission  coupled  to  the  cavity 
mode  goes  in  one  direction.  Using  equation  (14),  we  have: 


leril2=e2Gz 


which  gives  us: 


Pt=P2  =  fyErlp(E)y 


e2GL‘  - 1 


2G 


(50.a) 


and: 


P3=P4=ftrErjp(E)y 


e2GU- 1 


2G 


(50.b) 


Putting  (49)  and  (50)  together,  we  find: 


Pout=fyErjp(E) 


Yko 


[  K- 1  s  I  2(2Kb-K)  j 


(Cj+Cj) 


e2GL|  - 1 


2G 


+  (C3+C4) 


^-lll 


2G 


\ 


(51) 


Finally,  let  us  notice  that  this  is  Pout(E).  If  we  want  Po«(X),  since  E=hcA,  we  have  to  use: 


PouA)  =  PomCE)  I  1  =  Pout(E)^- 


The  factor  between  square  brackets  in  (51)  has  to  be  examined  carefully.  It  can  be 
shown  [17]  that  in  a  periodic  waveguide,  y  becomes  dependent  on  wavelength  and  con- 

m 

tains  a  factor  that  cancels  the  denominator.  The  expression  between  square  brackets  is 
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then  equal  to  the  y  in  an  uniform  waveguide  times  a  factor  of  order  unity  depending  on 
wavelength.  In  this  paper,  we  will  neglect  that  factor,  although  it  modifies  somewhat  the 
spectrum  [17]. 

Now,  only  two  factors  remain  to  be  specified  before  we  can  compute  the  emission 
spectrum.  These  are  the  spontaneous  emission  rate  r^(E)  and  the  GaAs  material  gain 
g,(E).  Both  will  be  functions  of  the  carrier  density.  For  the  spontaneous  emission  rate, 
we  know  [15]  that  for  the  non-thermal  equilibrium  case,  we  have: 


r»p(E)= 


_J _ 

trW 


-E2g(E,N) 

erE-A£PVk,T_  j 


(52) 


where  v  is  the  light’s  velocity  in  the  material,  E  is  the  photon  energy,  kB  is  Boltzmann’s 

constant  and  T  is  the  temperature.  The  difference  between  the  electron  and  hole  quasi- 

Fermi  levels  (Aep)  can  be  computed  from  the  carrier  density  N  by  using  the  Fermi-Dirac 

distribution  and  the  parabolic  approximation  for  the  band  extrema.  Using  this  and  exper- 

* 

imental  absorption  curves  from  [15],  we  developed  a  semi-empirical  model  for  the  gain 
as  a  function  of  E,  N,  T  and  NA,  where  NA  is  the  acceptor  concentration  (zinc  gives  a  p- 
type  doping).  Figure  17(a)  shows  typical  gain  curves  for  different  carrier  densities,  at 
constant  temperature,  and  figure  17(b)  for  different  temperatures,  at  constant  carrier  den¬ 
sity. 

At  this  point,  we  have  not  yet  incorporated  rate  equations  in  our  model  to  allow  us 
to  compute  the  L-I  characteristic.  Therefore,  we  will  postulate  a  carrier  density  and  com¬ 
pute  the  shape  of  the  emission  spectrum,  but  we  will  not  try  to  get  the  absolute  values. 
Thus  we  do  not  need  to  wony  about  the  constants  in  (50)  and  (51).  The  final  result  is 
shown  in  figure  18(a)  and  can  be  compared  with  the  experimental  spectra  shown  in  figure 
15  and  18(b).  The  resolution  of  the  measurements  in  figure  18(b)  is  0.7  nm.  Figure  18(b) 
was  measured  at  much  higher  currents,  in  pulsed  mode  operation.  We  see  that  the  theory 
and  the  experiment  are  similar  (with  reasonable  values  of  the  parameters),  but  that  the 
relative  intensities  of  the  peaks  are  not  reproduced.  This  may  be  due  to  the  fact  that  our 
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modei  is  linear  and  cannot  go  above  threshold  since,  as  mentioned  above,  we  do  not  use  a 
rate  equation.  Another  problem  is  that  we  assumed  a  uniform  carrier  density  in  the 
pumped  section.  This  is  most  probably  not  true  because  the  carrier  injection  is  not  uni¬ 
form  and  the  field  in  the  cavity  is  not  uniform  either.  The  field  is  maximum  at  the  phase 
shifter  and  decreases  exponentially  above  and  below  (see  [16]).  This  will  introduce  a 
non-uniformity  in  the  carrier  distribution  that  cannot  be  compensated  for  by  the  carrier 
diffusion  because  of  the  layered  structure.  So,  in  the  regions  of  low  field,  we  will  have 
more  spontaneous  emission.  Also,  the  bottom  of  the  mesa,  which  is  a  low  field  region,  is 
much  wider  than  the  top  (fig.  13(b)).  This  means  that  the  light  emitted  on  the  sides  sees  a 
different  vertical  structure  and  may  alter  the  emission  spectrum. 

Notice  in  figure  18  the  slight  shift  in  peak  position  with  increasing  carrier  density. 
This  is  due  to  the  change  in  the  index  of  refraction  of  the  GaAs  layers  with  the  carrier 
density.  In  our  model,  we  used  [6]: 

An  =  2.2xlCT®NI/3  -  1.6xlO-,4N2/3r  9-75x10“^  (20) 

where  An  is  the  index  change  and  N  is  the  carrier  density  in  cm-3.  We  also  incorporated  a 
dependence  on  temperature  [15]: 

An = 4x  1 0-4  (T-300)  (21) 

where  T  is  the  absolute  temperature  in  Kelvins.  We  assumed  that  the  change  in  index  of 
refraction  was  the  same  for  the  GaAs  and  the  GaAlAs  layers. 

This  effect  of  the  temperature  on  the  index  of  refraction  is  responsible  for  the  shift 
in  peak  positions  that  we  observe  in  figure  19,  which  shows  the  emission  spectrum  at 
several  temperatures  (with  the  carrier  density  kept  constant).  Figure  19(a)  shows  the 
theoretical  calculations  and  figure  19(b)  shows  the  experimental  result.  Once  again,  the 
agreement  is  good,  except  for  the  relative  peak  amplitudes.  We  notice  that  when  the 
temperature  changes,  the  relative  peak  amplitudes  change  also.  This  is  due  in  part  to  the 
modification  of  the  gain  spectrum  with  temperature,  as  shown  in  figure  17(b).  This 
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modificadon  is  due  to  a  number  of  factors,  the  main  ones  being  the  change  of  the 
bandgap  and  that  of  the  density  of  states  with  temperature.  We  see  that  the  wavelength 
of  maximum  gain  changes  and  so  does  the  value  of  that  maximum  gam.  When  the  mode 
position  coincides  with  the  maximum  gain,  the  peak  amplitude  is  maximum.  This  occurs 
around  300K  for  the  peak  at  884  nm  (at  300K).  When  we  change  the  temperature,  the 
maximum  gain  and  the  mode  shift  at  different  rates  and,  as  a  result,  the  peak  amplitude 
decreases. 

IX.  Conclusions 

In  this  paper,  we  have  developed  a  way  to  compute  the  propagation  constants  of  the 
modes  in  a  one-dimensional  periodic  structure  over  a  certain  band  of  frequency,  when 
material  dispersion  is  present  The  correct  interpretation  of  the  solutions  is  based  on  the 
direction  of  the  power  flow.  This  led  us  to  discover  that,  for  a  structure  with  even  sym¬ 
metry,  the  backward  and  forward  propagating  modes  exchange  their  roles  when  the  gain 
changes  sign  in  the  stop-band.  These  results  were  checked  in  the  two  layers  case  by 
obtaining  an  exact  solution.  We  then  derived  the  formulas  for  the  reflection  and 
transmission  coefficients  at  the  interface  between  a  periodic  and  a  uniform  section. 

Using  these  results,  we  developed  a  model  to  compute  the  emission  spectrum  of  the 
surface-emitting  laser  diode  (SELD)  below  threshold.  We  made  several  assumptions  in 
this  model.  The  main  ones  are:  (1)  The  Bloch  wave  description  for  the  electric  field  is 
accurate.  (2)  The  carrier  density  is  uniform  in  the  pumped  region.  (3)  The  horizontal 
cross-section  (parallel  to  the  layers)  is  constant  from  the  top  to  the  bottom.  We  think  that 
the  first  assumption  is  a  good  approximation,  but  that  the  second  and  third  ones  are  par¬ 
tially  responsible  for  the  discrepancies  between  the  theory  and  the  experimental  results. 
Our  confidence  in  the  first  assumption  comes  from  the  good  agreement  between  theory 
and  experiment  for  the  reflectivity  of  the  layered  structure  (section  VI). 
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As  mentioned  in  section  VIII  and  in  the  previous  paragraph,  the  assumption  of  a 
uniform  carrier  density  is  probably  not  a  good  one.  A  variation  in  carrier  density  will 
result  in  a  longitudinal  gain  non-uniformity,  which  may  change  the  threshold  gain  for  the 
modes  and  hence  explain  the  difference  in  the  relative  peak  amplitudes  between  theory 
and  experiment.  The  fact  that  the  cavity  is  really  trapezoidal  (and  not  rectangular,  as 
implied  by  the  third  assumption)  may  give  a  significant  contribution  to  spontaneous 
emission  and  modify  the  spectrum,  since  the  light  emitted  on  the  sides  sees  a  different 
vertical  structure  (see  fig.  13). 

In  conclusion,  we  can  say  that  our  model  is  sufficient  to  explain  a  lot  of  the 
observed  characteristics  of  the  spectrum.  To  get  a  better  agreement,  we  need  to  extend 
its  validity  above  threshold.  To  do  that,  we  need  to  incorporate  rate  equations  for  the  car¬ 
rier  density.  The  next  step  would  be  to  take  into  account  the  longitudinal  non-uniformity 
of  the  carrier  density,  but  that  is  much  more  difficult  and  can  probably  not  be  done  using 
the  theory  developed  in  this  paper.  One  way  to  do  it  would  be  to  use  a  transfer  matrix  for 
each  period  (see  for  example  [13])  and  solve  self-consistcntly. 

On  the  experimental  side,  we  are  now  working  on  several  improvements,  such  as 
using  reactive  ion-beam  etching  to  get  vertical  side-walls  to  the  cavity,  which  would 
avoid  the  problems  due  to  a  trapezoidal  cavity.  We  are  considering  replacing  the  zinc 
diffusion,  which  is  difficult  to  control,  by  a  two  steps  LPE  process.  We  grow  first  n-type 
GaALAs  on  one  side  of  the  mesa,  then  p-type  GaAlAs  on  the  other  side,  so  that  we  still 
get  a  transverse  p-n  junction.  We  are  also  working  on  distributed  Bragg  reflector 
(DBR)-type  structures,  where  we  have  a  relatively  long  uniform  active  region  (3-5  pm) 
between  a  bottom  unpumped  multilayer  and  a  top  dielectric  multilayer. 
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Appendix  A 


In  this  appendix,  we  derive  the  exact  expression  of  the  average  over  one  period  of 
the  rime  averaged  Poynting  vector.  We  start  with  an  electric  field  of  the  form: 

Ey=  [ e>m  (A.l) 

Using  Maxwell’s  equations,  this  gives: 

Hx= — —(g-jP)f  e***  (A. 2) 

O&Po  l 

The  time  averaged  Poynting  vector  is: 

S  -  jRe(gxi?) = -  jRe(EyH;>^  (A.3) 

Which  gives  us: 

Sz=  -P  -  I Cj 1 2 e2**—  lc2l2e“2*z  +— 2—  im  (A.4) 

20)Po  l  J  a>Po  L 

Now,  for  Ey=A(z)erz,  where  A(z)  is  given  by  (27),  the  average  of  Sz  over  one  spatial 
period  is  given  by: 

1  * 

<s*>=—  l Szdz 


Pi  i  -e-2*’*  ,  ,  2  e2*’1*  - 1 

T  ~li^+ISl'  “5 T 


(A.5) 
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Appendix  B 

This  appendix  gives  the  expressions  for  the  reflection  and  transmission  coefficients 
anc  for  the  propagation  factors  used  to  compute  the  reflectivity,  in  (44),  and  the  emission 
spectrum,  in  (46>-(51).  With  Li=z^~z,  and  L3=  zV-z'j  ivfig.  1 1  and  16),  we  have: 
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(B.2) 


The  coefficients  sf  and  s*,  are  defined  in  (34).  The  numerical  index  indicates  at  which 
position  Zj  they  are  computed  (and  ko=2jt/X). 

For  the  reflectivity  calculation  (44),  we  have  (fig.  11): 

1+Sb4 


R  j - Sf4 


1+Sf4 


Tj  =  (l  -Sf4sM) 
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1+Sf4 


(B.3) 


Notice  that  for  the  emission  spectrum  calculation,  Lj  is  different  from  the  Lj  in  the 


reflectivity  calculation.  For  the  spectrum  (46)-(51),  we  get  (fig.  16): 
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R5  =  R$  + 


t6t7p5P6R8 

i-r7r8p5p6 


(B.4) 


where  R5  is  an  effective  reflection  coefficient,  for  everything  to  the  left  of  the  L3/L4  inter¬ 
face  in  figure  16.  And  we  use  (with  L4  =  z's~z'4): 
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The  superscript  u  for  s  indicates  here  that  it  is  the  s  of  the  unpumped  section  (L4)  that  has 
to  be  used. 
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Figure  captions 

Figure  1:  Examples  of  periodic  structure  that  give  a  one-dimensional  wave  equation:  (a) 
waveguides,  (b)  multilayers. 

Figure  2:  Definition  of  the  parameters  for  the  two-layers  structure. 

Figure  3:  Intensity  gain  in  GaAs,  as  a  function  of  free  space  wavelength  for  three  dif¬ 
ferent  carrier  densities:  (1)  N=0.2xl018  cm-3,  (2)  N=1.0xl018  cm-3,  (3) 

N=2.0xiOiS  cm-3. 

Figure  4:  G  and  5^  for  the  structure  of  table  1:  (a)  without  gain  or  loss,  (b)  with  2000 
cm-1  average  gain,  but  neglecting  gain  periodicity,  (c)  as  in  b,  but  with  gain  periodi¬ 
city  (see  text). 

Figure  5:  G  and  8^  for  the  structure  of  table  1  and  the  intensity  gains  of  figure  3, 
corresponding  to  three  different  carrier  densities:  (a)  N=0.2xi018  cm-3,  (b) 

N=1.0xi018  cm_3,Xc)  N=2.0xl018  cm-3. 

Figure  6:  is)  as  a  function  of  wavelength  for  the  two  solutions  of  eq.  (13),  at  three  diffrem 
carrier  densities  (Fig.  3):  (a)  N=0.2xl018  cm-3,  (b)  N=1.0xl018  cm-3,  (c) 

N=2.0xi018  cm-3. 

Figure  7:  Re(s)  and  Im(s)  as  functions  of  wavelength,  at  three  different  carrier  densities 
(Fig.  3):  (a)  N=0.2xl018  cm-3,  (b)  N=1.0xl018  cm"3,  (c)  N=2.0xl018  cm"3. 

Figure  8:  Definition  of  the  transmission  and  reflection  coefficients  at  a  periodic/unifonn 
regions  interface  (equation  (39)). 

Figure  9:  Definition  of  the  transmission  and  reflection  coefficients  at  a  phase  shifter 
(equation  (43)). 

Figure  10:  Model  used  for  reflectivity  and  emission  spectrum  calculations.  One  pair  con¬ 
sists  of  a  quater-wavelength  layer  of  GaAs  and  a  quarter-wavelength  layer  of  GaA- 
LAs.  Dimensions  are  given  in  Table  1. 
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Figure  1 1:  Definition  of  the  coefficients  for  the  reflectivity  calculation  (44). 

Figure  12:  (a)  Result  of  the  reflectivity  calculation,  with  a  center  wavelength  of  Xq=  893 
nm.  (b)  Comparaison  of  experiment  with  best  theoretical  result  (see  text). 

Figure  13:  (a)  Structure  of  the  surface-emitting  laser  diode  (SELD).  (b)  Scanning  elec¬ 
tron  microscope  picture  of  a  SELD. 

Figure  14:  L-I  characteristic  of  a  SELD  at  room  temperature,  CW  operation. 

Figure  15:  Emission  spectrum  of  a  SELD  at  different  drive  currents,  room  temperature, 
CW  operation. 

Figure  16:  Definition  of  the  coefficients  for  the  emission  spectrum  calculation  (equations 
(46-51). 

Figure  17:  GaAs  intensity  gain  curves  (NA  =  2xl0ls  cm*3):  (a)  at  T=300K,  for  several  car¬ 
rier  densities,  (b)  at  N=l.Qxl018  cm-3  and  several  temperatures.  These  were 
obtained  with  the  same  model  as  for  figure  3. 

Figure  18:  Emission  spectrum  for  several  carrier  densities  (currents)  at  300K:  (a)  theory 
(Xq=  890  nm,  phase  shifter  thickness  =  0.4  A),  (b)  experiment  (pulsed  mode  opera¬ 
tion). 

Figure  19:  Emission  spectrum  for  several  temperature,  at  constant  carrier  density 
(current):  (a)  theory  (N=  l.OxlO1*  cm-3,  same  parameters  as  in  figure  18),  (b)  experi¬ 
ment  (1=320  mA,  pulsed  operation). 
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On  The  Emission  Properties  of  Surface-Emitting  DFB  and  DBR 
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We  use  a  single  mode  transfer  equation  model  to  study  the 
performance  of  vertical-cavity  surface-emitting  DFB  and  DBR 
semiconductor  lasers  above  threshold.  We  find  that  DBR  lasers  exhibit  less 
envelope  spatial  hole  burning  than  both  single  and  double  phase-shifted 
DFB  lasers  and  therefore  may  achieve  more  stable  single  longitudinal  mode 
operation. 


Surface-emitting  semiconductor  lasers  are  being  developed  for  two 
dimensional  array  applications  and  may  eventually  be  useful  in 
optoelectronic  integrated  circuit  technology.  The  designs  can  be 
distinguished  according  to  whether  the  resonator  feedback  is  in  the  plane 
of  the  wafer* or  normal  to  it~*'^;  Fabry-Perot,  distributed  feedback 
(DFB),  distributed  Bragg  reflector  (DBR),  and  mixed  Fabry-Perot-DBR-type 
devices  have  been  demonstrated.  The  threshold  current  of  these  lasers 
must  be  minimized  to  reduce  heating  effects,  total  array  power 
requirements,  and  the  size  of  on-chip  drive  electronics,  but  a  high  output 
efficiency  is  also  desirable.  A  model  intended  to  determine  the  optimum 
structure  should  be  capable  of  calculating  the  output  characteristics  above 
threshold.  Several  authors  have  used  coupled  mode  theory  to  analyze  DFB 
lasers  in  the  above  threshold  regime.  Hill  and  Watanabe^  constructed  a 
numerical  treatment  and  considered  the  problems  caused  by  parasitic 
losses.  Haus*®  derived  analytical  formulas  for  the  output  power  of 
overcoupled  lasers,  and  Szczepanski *  *  extended  these  results  to  the 
undercoupled  case.  These  authors  did  not  consider  the  effect  of  phase  shift 
layers,  and  the  effect  of  facet  reflections  was  considered  only  in  Ref.  11. 
Also  these  results  are  more  appropriate  to  edge-emitting  lasers;  the 
coupled  mode  formalism  assumes  weak  coupling  between  the 
counterpropagating  waves,  and  surface  emitting  lasers  nominally  operate 
in  the  strong  coupling  regime.  Typically  the  distributed  feedback 
coefficient  K  ~  5000-10000  cm‘^  for  surface  emitters,  and  K  ~  50  cm'*  for 
edge  emitters.  In  order  to  take  this  strong  coupling  effect  into 
consideration  and  to  include  the  effects  of  gain  saturation,  we  have 
extended  the  discrete  propagation  techniques  of  Bjork  and  Nilsson*  ^  and 
have  applied  this  model  to  a  comparison  of  vertical  cavity  surface-emitting 


DFB  and  distributed  Bragg  reflector  (DBR)  lasers.  We  also  analyze  the 

n 

performance  of  a  structure  we  have  recently  developed  .  The  methods  of 
Ref.  12  are  particularly  amenable  to  computer  calculation  and  well  suited 
to  analyzing  resonators  having  a  multilayer  structure.  They  are  also 
flexible  enough  to  be  used  in  cases  involving  facet  reflections,  phase  shift 
layers,  amplifying  regions,  and  dielectric  coalings. 

The  type  of  laser  we  will  consider  employs  a  multilayer 
heterostructure  for  distributed  feedback  and  a  buried  heterostructure 
waveguide  for  optical  mode  confinement.  Phase  shift  layers  or  amplifying 
layers  may  be  placed  anywhere  along  the  resonator  length.  The  equations 
we  use  to  model  such  a  laser  comprise  two  sets.  The  first  set  treats  the 
free  propagation  of  the  forward  and  backward  waves  within  a  layer 
according  to 

AjCZj+lj)  =  AjCzj)  expfiPjij  +  (1/2)  Qj(Aj,  BjJljJ  (1) 

BjUj+lj)  =  Bj(Zj)  expHPjlj  -  (1/2)  gi(Aj,  B;)l,l  (2) 

9 i< Aj,  Bj)  .  g0/(1  +  |Aj(Zj)|2  +  IBj(Zj)l^)  (3) 

where  A-  and  Bj  are  the  normalized  forward  and  backward  wave 
amplitudes  in  the  i-th  layer,  Zj  is  the  axial  position  of  layer  i,  Pj  is  the 
propagation  constant,  lj  is  the  i-th  layer  thickness,  and  gj  and  gQ  are  the 

saturated  and  unsaturated  intensity  gain  coefficients  respectively.  In  Eq. 
(3)  we  have  assumed  that  carrier  diffusion  effectively  eliminates  any 
longitudinal  spatial  hole  burning  within  a  layer;  we  also  ignore  the 
transverse  variation  in  the  laser  mode  and  homogeneously  saturated  gain. 
The  fields  are  propagated  across  an  interface  between  layers  i  and  j 


(Pi -Pi)  (Pj  +  Pj) 

Bi  -  Fii^— + 


p«  • 

Hjj  =  Z  J  Jdxdy  WjWj 


(5) 

(6) 

(7) 


where  Wj  is  the  waveguide  mode  field  function  in  layer  i,  and  the  sum  in 

Eq.  (7)  is  over  polarization.  Eqs.  (4-7)  have  been  derived  following  results 
by  Wang1-*.  The  overlap  integrals  in  Eq.  (6)  account  for  the  mismatch  in 
mode  profile  between  two  layers.  It  is  assumed  that  the  power  scattered 
into  higher  order  waveguide  modes  is  negligible;  this  is  a  good 
approximation  for  the  surface  emitting  lasers  demonstrated  to  date.  In  the 
case  of  an  interface  between  a  waveguiding  region  and  a  non-waveguiding 
region,  we  assume  that  to  first  order  and  for  short  propagation  distances 
the  transverse  mode  profile  in  the  free  space  region  is  the  same  as  that  in 
the  waveguide;  we  also  use  the  plane  wave  propagation  constant  for  the 
non-waveguide  region.  Eqs.  (1*7)  are  used  repetitively  to  calculate  the 

field  amplitude  in  all  layers;  the  boundary  conditions  are  satisfied  only  for 
a  particular  set  of  values  of  g@,  X,  Bj,  and  An+2  (n  is  the  number  of  layers). 

We  use  a  quasi-Newton  optimization  routine  to  efficiently  calculate  the 
solution  for  each  value  of  the  output  power1 

We  show  in  Fig.  1  the  forward  and  backward  output  power  versus 
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pumping  calculated  for  the  lowest  order  longitudinal  mode  (at  882.6  nm) 

of  the  structure  of  Ref.  7.  The  resonator  consists  of  80  pairs  of 

A I  q  ^Guq  yAs/GaAs  tuned  to  883  nm.  The  quarter-wave  phase  shift  layer 

is  used  to  ensure  oscillation  at  the  Bragg  mode;  its  placement  20  pairs 
from  the  output  facet  effectively  diverts  most  of  the  power  in  that 
direction.  The  slight  superlinearity  in  the  output  curve  can  be  explained 
by  examining  the  axial  distribution  of  the  intensity,  plotted  in  Fig.  2a,  at  a 
pumping  twice  the  threshold  pumping.  It  is  seen  that  the  intensity  is 

maximized  at  the  phase  shift  layer,  and  the  gain  is  saturated  only  in  the 

front  portion  of  the  resonator.  This  mode  will  not  tend  to  oscillate 
coherently  since  spontaneous  emission  processes  dominate  in  the  back 

portion  of  the  resonator.  In  fact  the  gain  there  is  available  for  extraction 
by  the  next  mode  with  shorter  wavelength  (at  860.8  nm),  since  its 
intensity  is  maximized  in  that  section  of  the  structure  (Fig.  2b).  For  this 

laser  the  threshold  gain  of  the  mode  at  860.8  nm  is  roughly  four  times  that 

of  the  Bragg  mode,  and  therefore  its  contribution  to  the  output  spectrum 
will  likely  be  as  a  broad  amplified  spontaneous  emission  background.  A 
similar  effect  was  noted  by  Usami  et  al.16  in  a  study  of  edge-emitting 
lasers;  the  large  coupling  coefficients  used  here  (K  ~  5400  cm'1,  KL  ~  5.4) 
exacerbate  this  problem.  The  large  oscillation  linewidth  and  multimode 
spectrum  observed  in  Ref.  7  may  possibly  be  explained  by  these 
considerations. 

It  has  been  suggested  recently16'18  that  the  use  of  multiple  phase 

shift  layers  can  help  reduce  the  problems  associated  with  spatial  hole 

burning.  We  examine  by  way  of  example  a  distributed  feedback  design 
with  75  pairs  of  A1q  ^GaQ  j  As/GaAs  tuned  to  839  nm,  and  with  phase  shift 

layers  inserted  at  the  25-th  and  55-th  pairs  from  the  back  facet.  Both 


shifters  are  0.425  waves  in  thickness,  and  the  estimated  KL  for  this 
structure  is  12.5.  These  parameters  were  chosen  to  smooth  the  saturated 
gain  distribution  as  well  as  force  lasing  near  the  GaAs  spectral  gain 
maximum.  The  axial  intensity  distribution  for  the  first  lasing  mode  at  875 
nm  is  shown  in  Fig.  3a,  again  at  twice  the  threshold.  The  intensity  is  a 
maximum  at  both  phase  shift  layers,  and  the  fraction  of  the  resonator 
experiencing  gain  saturation  and  coherent  emission  is  larger  than  in  the 
single  phase  shifter  design.  The  large  detuning  of  the  mode  from  the  Bragg 
wavelength  aids  in  smoothing  the  intensity  distribution  between  the  phase 
shift  layers.  Unfortunately  this  mode  is  still  susceptible  to  gain 
competition  from  neighboring  modes  as  is  confirmed  by  examining  the 
axial  distribution  of  the  higher  threshold  mode  at  892  nm  (Fig.  3b).  This 
mode  can  extract  the  spatial  gain  that  is  unused  by  the  main  lasing  mode. 

Simple  DBR  lasers  do  not  in  general  suffer  from  these  problems.  The 
boundary  conditions  imposed  by  the  two  endmirrors  cause  the  axial 
intensity  to  be  a  maximum  in  the  amplifying  layer,  leading  to  efficient  gain 
conversion.  This  is  illustrated  in  Fig.  4  for  a  representative  DBR  structure 
(for  one  endmirror  KL~3.50,  for  the  other  KL~1.66;  the  amplifying  layer  is 
2  pm  thick).  For  the  higher  order  longitudinal  modes  the  increased 
detuning  from  the  Bragg  wavelength  causes  the  intensity  distribution  to  be 
a  maximum  in  the  DBR  reflector  and  outside  the  amplifying  layer,  reducing 
the  gain  extraction  and  ensuring  good  mode  discrimination.  Mixed  DBR- 
Fabry-Perot  lasers  have  exhibited  high  mode  selectivity  and  narrow 
linewidth^’^,  lending  credence  to  these  conclusions.  Our  theory  can  be 
used  to  model  these  lasers  by  assuming  plane  wave  propagation  in  the 
dielectric  mirror,  as  discussed  above.  A  prediction  of  the  oscillation 
linewidth  of  single  mode  emitters  is  outside  the  scope  of  our  present 


treatment,  as  it  is  determined  in  part  by  factors  such  as  layer  thickness 
fluctuations  which  we  have  not  included. 

In  conclusion,  we  have  applied  discrete  propagation  methods  to 
analyze  the  emission  properties  of  DFB  and  DBR  surface  emitting  lasers. 
For  DFB  lasers  with  either  one  or  two  phase  shifting  layers  the  mode 
discrimination  can  be  rather  poor  due  to  competition  for  the  spatially 
varying  gain.  Our  results  suggest  that  DBR  lasers  are  less  susceptible  to 
these  spatial  hole  burning  effects. 

This  research  was  supported  in  part  by  the  Office  of  Naval  Research 
and  the  Air  Force  Rome  Air  Development  Center. 
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Figure  Captions  \ 

Fig.  1.  Power  output  versus  pumping  parameter. 

Fig.  2.  Single  phase-shift  laser  axial  intensity  distribution  for  the  Bragg 

mode  (a),  and  next  higher  order  longitudinal  mode  (b)  . 

Fig.  3.  Double  phase-shift  laser  axial  intensity  distribution  for  the  lowest 

order  (a)  and  next  to  lowest  order  (b)  axial  modes. 

Fig.  4.  Axial  intensity  distribution  of  the  fundamental  mode  of  a  DBR  laser. 
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ABSTRACT 

The  structure  and  performance  characteristics  of  a  GaAs  surface  emitting 
laser  (SEL)  diode  and  a  2-D  array  of  SELs  is  reported.  The  SEL  has  a  very 
short  optical  cavity  with  a  distributed  feedback  structure,  and  emit3  light 
perpendicular  to  the  wafer  surface.  The  laser  shows  low  threshold  current  and 
has  an  estimated  high  frequency  response  of  tens  of  GHz.  A  2-D  array  of  SELs 
with  6  x  16  elements  was  designed  and  demonstrated. 


1-.  INTRODUCTION 

A  two  dimensional (2-D)  array  of  laser  diodes  i3  desirable  for  many 
applications  including  optical  signal  processing,  optical  interconnect,  and 
optical  computing.  In  particular,  a  2-D  array  with  high  modulation  rate 
capability  and  individually  addressable  control  can  have  application  in 
microwave  optics  3uch  as  an  optically  controlled  phased  array  microwave 
system.  In  this  paper  we  report  on  a  new  type  of  2-D  array  made  of  surface 
emitting  laser  diodes  we  demonstrated  earlier.1  The  SEL  diode  has  many  unique 
features  including  vertical  emission,  short  optical  cavity,  low  threshold 
current,  and  the  potential  for  very  high  rate  modulation.  Mo3t  importantly,  it 
can  be  scaled  into  a  2-D  array  easily  without  the  rack  and  stack  technique 
used  in  assembling  conventional  laser  diode  arrays. 

In  the  following  sections,  we  briefly  describe  the  structure  of  this  laser 
diode  and  present  the  key  performance  characteristics.  The  modulation  response 
of  the  SEL  diode  is  also  discussed.  Finally  we  illustrate  the  configuration 
and  demonstration  of  a  2-D  laser  array  made  of  the  SEL  diodes. 


2.  LASER  DIODE  STRUCTURE 


A  schematic  sketch  of  the  SEL  is  shown  in  Figure  1.  The  laser  diode 
structure  is  a  surface  emitting  type  in  which  the  light  emission  is 
perpendicular  to  the  wafer  surface  instead  of  parallel  to  the  wafer  surface  as 
in  conventional  edge  emitting  laser  diodes.  The  SEL  utilizes  a  transverse 
buried  hetero junction  structure  to  inject  current  laterally  and  to  confine  the 
injected  carriers,  and  has  a  distributed  feedback  (DFB)  vertical  cavity  type 
resonator.  The  active  region  is  a  multilayer  structure  which  is  made  of 
multiple  pairs  of  Alo.3Gao.7As  and  GaAs  layers  with  quarter-wavelength  optical 
thickness  per  layer.  Typically,  the  multilayer  structure  has  20  pairs  of  these 


thin  layers  at  the  top  and  60  pairs  at  the  bottom,  and  a  GaAs  phase  shifting 
layer  of  about  quarter-wavelength  thickness  in  between.  The  phase-  shifting 
layer  allows  the  laser  to  lase  at  the  Bragg  wavelength  of  the  multilayer 
structure.  This  multilayer  structure  also  acts  as  a  DFB  optical  cavity  for  the 
laser  diode.  The  buried  heterostructure  with  a  transverse  junction  is  formed 
by  the  N-  and  ?-AlGaAs  hetero junction  cladding  layers  surrounding  the  active 
region . 

Because  of  the  lower  turn-on  voltage  for  the  GaAs  junction  than  that  of 
the  AlGaAs  junction,  the  carriers  are  predominantly  injected  into  the  GaAs 
part  of  the  active  region.  The  radiative  recombination  of  these  carriers  in 
the  active  region  and  the  feedback  effect  of  the  DFB  cavity  produce  vertical 
laser  emission  perpendicular  to  the  junction  and  the  wafer  surface  as  shown  in 
Fig.  1. 


Figure  1.  Schematic  diagram  of  surface  emitting  laser  diode 


3.  TYPICAL  SEL  PERFORMANCE  CHARACTERISTICS 


A  typical  SEL  ha3  a  rectangular  emitting  area  of  about  3  x  15  um  and  an 
electrode  on  each  side  (p-  and  n-junction)  providing  for  lateral  current 
injection.  With  this  type  of  structure  and  arrangement,  the  SEL  can  be  tested 
as  grown  on  the  wafer  without  the  cleavage  process  required  for  the 
conventional  laser  diodes.  It  also  makes  2-D  array  formation  and 
optoelectronics  integration  easier.  The  following  performance  characteristics 


of  SEL  were  measured  on  the  wafer  basis  and  under  the  normal  room-temperature 
operation . 


One  of  the  distinct  cha  racter  ist  ics  of  this  SEL  is  its  emission  pattern. 
Figure  2  shows  a  near-field  emission  pattern  as  viewed  with  a  photomitr -'scope . 
The  emission  pattern  is  tightly  confined  within  the  rectangular  emi t t ing 
surface  area,  and  is  fairly  uniform  across  the  area  as  shown  by  the  intensity 
profile.  The  far-field  typically  shows  a  nearly  circular  beam  pattern  with  a 
beam  divergence  angle  of  about  S  deg.  The  laser  emission  wavelength  of  this 
SEL  is  typically  around  884  nm  near  the  GaAs  bandgap . 


Figure  2.  Near-field  emission  pattern 

Another  characteristic  of  this  SEL  that  is  very  important  for  high  speed 
modulation  and  optical  signal  processing  is  the  relatively  low  threshold 
current  operation  requirement  and  linear  light  output  dependence  on  the 
driving  current.  Figure  3  shows  a  typical  light  output  versus  driving  current 
relation  of  the  SEL  tested  under  cw  operating  conditions.  The  light  output 
increases  linearly  with  the  current  up  to  25  mA.  This  linearity  is  desirable 
for  distortion  free  analog  modulation.  The  output  power  of  the  SEL  is  about 
90  uW  at  25  mA  for  this  case,  a  higher  output  of  up  to  0.5  mW  was  obtained  at 
higher  operating  current.  The  threshold  current  is  about  2  mA  based  on  this 
data,  which  is  substantially  lower  than  the  conventional  laser  diode.  This  low 
threshold  current  characteristic  together  with  the  short  laser  cavity 
configuration  promise  a  very  high  speed  modulation  capability. 


I 


0 


5 


i  i  > - 1  •  i - ' - 1 

10  15  20  25 

DC  CURRENT  (mA) 


Figure  3.  Light  output  vs  driving  current  characteristics 


4.  MODULATION  RESPONSE  QF  SEL 


The  modulation  bandwidth  of  a  laser  diode  is  typically  limited  by  the 
relaxation  oscillation  of  the  laser  diode.  The  relaxation  oscillation 
frequency,  f,  i3  related  to  the  laser  device  parameters  through  the  following 
approximate  expression: 


Where  Te  and  Tp  are  lifetime  of  carriers  and  photons,  Jo  is  the  threshold 
current  density  and  J  is  the  operating  current  density.  For  the  SEL,  the 
threshold  current  is  about  2  mA  and  the  cavity  length  is  about  5  pm.  If  we 
assume  the  confinement  factor  is  0.5  ,  internal  loss  is  10  cm-1,  a  cavity  of 
effective  length  5  nm,  and  reflectivity  of  70  %,  then  the  photon  lifetime  i3 
about  8  x  10~14  sec.  For  a  carrier  lifetime  of  2  nsec  (  typical  for  GaAs  ), 
the  relaxation  frequency  is  about  25  GHz  for  an  operating  current  of  10  mA. 
Therefore,  high  speed  modulation  can  be  expected  from  this  type  of  SEL. 


Preliminary  investigation  of  the  pulse  response  characteristics  of  the  SEL 
diode  was  conducted,  using  a  50  nsec  current  pulse  at  1  KHz  repetition  rate 
and  a  40  psec  rise  time  pnotcdetector .  Figure  4  shows  the  initial  result  of 
the  light  output  response  to  the  current  pulse  taken  by  a  sampling  scope  of 
20  GHz  bandwidth.  From  this  initial  data,  the  rise  time  of  the  light  emission 
seems  to  follow  that  of  the  current  pulse,  which  has  a  risetime  of  about  5 
nsec.  A  slight  relaxation  oscillation  is  observable  on  the  light  output 
response  but  is  not  very  pronounced.  Further  investigation  using  faster  rise 
time  current  pulse  is  planned.  In  addition,  a  new  improved  SEL  structure  to 
reduce  the  parasitic  capacitance  and  series  resistance  is  contemplated. 
Further  results  on  this  will  be  reported  later. 
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Figure  4.  Response  of  light  output  to  current  pulse 


5.  SEL  ARRAY  CONFIGURATION  AND  EMISSION 

Because  of  the  planar  structure  and  surface  emission  of  SEL,  formation  of 
1-D  and  2 -D  array  is  relatively  simple  to  accomplish.  For  example,  a  3  x  3  of 
SEL  array  can  be  formed  easily  at  illustrated  schematically  in  Figure  5.  By 
properly  design  of  the  masks  for  the  array,  the  SEL  array  can  be  fabricated 
using  the  same  teguence  of  processing  step  as  the  fabrication  of  individual 
SELs  . 
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Figure  5.  Schematic  configuration  of  a  3  x  3  SEL  array 


Using  the  DFB  type  SEL  diode  structure  described  previously,  we  fabricated 
a  number  of  arrays  with  1-D  and  2-D  configurations.  Figure  6  shows  the 
photomicroscopic  picture  of  a  2-D  SEL  array  with  6  x  16  elements  on  a  GaAs 
chip  of  about  2  mm  x  2  mm.  The  light  emitting  area  is  adjacent  to  the  bright 
shinning  rectangular  region  of  Au  electrode  pads.  The  spacing  between  each  SEL 
element  is  about  100  pm,  limited  by  the  electrode  pad  size  which  is  about  80pm 
x  80  pm.  The  light  emitting  area  varies  from  3  pm  x  15  pm  to  3  pm  x  20  pm. 

The  array  chip  was  mounted  on  a  standard  multiple  pin  chip  carrier  for 
testing.  Figure  7  shows  the  emission  pattern  of  the  SEL  array  under  cw 
operation.  Because  of  the  different  emission  area  of  SEL  elements  in  the  array 
and  possible  nonuniformity  in  the  wafer  and  fabrication  process,  the  threshold 
current  for  each  SEL  varies  widely  between  a  few  mA  to  a  few  tens  of  mA. 
Therefore,  the  array  emission  pattern  was  taken  by  adjusting  the  driving 
current  of  each  SEL  to  provide  a  nearly  uniform  emission  intensity.  As  can  be 
seen  from  the  emission  pattern,  the  overall  array  has  a  fairly  good  yield  with 
about  95  %  of  the  SELs  in  the  array  emitting  radiation. 

This  array  configuration  allows  each  SEL  element  in  the  array  to  be 
addressed  individually  as  demonstrated  by  the  emission  pattern  shown  in  Figure 
7  and  provides  possibility  of  modulating  each  SEL  element  independently. 
Furthermore,  2-D  optoelectronic  integration  of  the  SEL  with  other  electronic 
components  such  as  a  MESFET  driver  to  provide  very  high  speed  operation  is 
al3o  possible  for  this  type  of  SEL. 
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Figure  6.  Photomicroscopic  view  of  a  6  x  16  SEL  array 
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We  report  a  new  type  of  2-D  laser  array  made  of  SEL  diodes.  The  SEL  diode 
has  a  transverse  buried  hetero junction  structure  with  a  short  distributed 
feedback  laser  cavity,  and  emits  vertical  radiation  perpendicular  to  the  wafer 
surface.  The  laser  shows  very  low  threshold  current  in  the  order  of  a  few  mA. 
These  key  features  together  with  the  preliminary  data  and  estimation  suggest 
that  the  SEL  and  its  array  can  be  modulated  at  a  very  high  rate  in  the  order 
of  tens  of  GHz.  We  describe  the  design  of  2-D  laser  array  based  on  the  SEL  and 
demonstrate  the  operation  of  a  6  x  16  2-D  SEL  laser  array.  The  result 
indicates  that  a  2-D  SEL  array  is  scalable  and  individually  addressable.  With 
many  unique  features  of  the  SEL  and  2-D  array,  we  believe  the  SEL  and  SEL 
array  will  be  very  attractive  for  use  in  microwave  optics,  optical  signal 
processing,  and  optical  interconnect  applications  in  the  future. 
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